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ACCELERATION OF SERIES 

by 
R.W. Gosper 

ABSTRACT 

The rate of convergence of infinite series can be accelerated by a suitable 
splitting of each term into two parts and then combining the second, part 
of the n-th term with the first part of the (n+l)-th term to get a new series 
and leaving the first part of the first term as an "orphan". Repeating this 
process an infinite number of times, the series will often approach zero, and 
we obtain the series of orphans, which may converge faster than the original 
series. Heuristics for determining the splits are given. Various mathematical 
constants, originally defined as series having a term ratio which approaches 
1, are accelerated into series having a term ratio less than 1. This is done 
with the constants of Euler and Catalan. The series for pi/4 = arctan 1 is 
transformed into a variety of series, among which is one having a term ratio 
of 1/27 and another having a term ratio of 54/3125. A series for 1/pi is 
found. which has a term ratio of 1/64 and each term of which is an integer 
divided by a power of 2, thus making it easy to evaluate the sum in binary 
arithmetic. We express zeta(3) in terms of pi+3 and a series having a term 
ratio of 1/16. Various hypergeometric function identities are found, as well 
as a series for (arcsin y)+2 curiously related to a series for y arcsin y. 
Convergence can also be accelerated for finite suras, as is shown for. the har- 
monic numbers. The sum of the reciprocals of the Fibonacci numbers has been 
expressed as a series having the convergence rate of a theta function. Finally, 
it is shown that a series whose n-th term ratio is (n+p) (n+q)/(n+r) (n+s) , where 
p, q, r, s are integers, is equal to c +• d pi+2, where c and d are rational. 
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Introduct i on 

This paper describes a way to manipulate sums to 
produce new ones which converge faster. For example, knowing- 
only that 

pi/4 - 1 - 1/3 + 1/5 - 1/7 + ... 
we can find that 

Pi 2 112 .13 4 15 6 

— - — (G {13 (20 (27 - , 

4 . 15 3 7 11 3 13 17 3 19 23 

without even using calculus.' Now the first series is a lot 
simpler, but it would require bi I I ions of terms to get 10 digi ts. 
of pi, whereas the second series only needs 7 terms, since 
each term is less than l/27th of the previous one. 
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Another example is the sum of the reciprocals of the cubes; 



1 1, 1 
zeta 3 » — + — + — + ... 

3 3 3 

12 3 n>l n 

Ue will show that this slow series 



Ev 



5 1 111 221 33 

- ( - - - - ( - - - . ( _ _ --(... 

A 1 342 563 78 



k-1 
5 V^ <-) 



E 



2 *— ' 3 2k 
k>l k ( ) 

where each term is less than 1/4 of the previous. 

The treatment will be informal with few, if any, proofs. 

The method consists of splitting each series term into two 
pieces, then recombining the pieces of adjacent terms to 
form a new series with smaller terms, and a 'leftover quanti ty 
called an "orphan". This process is repeated until a formula 
for the orphans is deduced, wherupon the shrinking series is 
discarded in favor of the series of orphans. The simplest 
case of this is an old trick known as 



Euler's Transfor 



m 



Although it is usually applied to a I ternat ing ser ies, Euler's 
transform is more neatly derived for the general series 



E 

n>8 



a ■a+a+a + a+... 
n . 8 1 2 3 



a a + a a + a a + a 

8 8 1 12 2 3 

- — + + + + 

2 2 2 2 



E 



a 

+ \ 'U a 
2 Z-/ n 

n>0 
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a + a 
n n+1 . 
where U is the averaging operator, U a « 

n 2 

Thus we get the orphan term (a /2) and a new series upon which 

2 

we may repeat this operation to get 



a U a ^ — > 2 
+ + 



2>n 



2 2 . n>0 

a + 2 a + a 
2 n n+1 n+2 . 
where U a - U (U a ) - : 

■ n . n ' '4 

If we do this k times we wi I I get 

. k-1 
a U a U a 

8 . ^-> k 

2 2 2 t-^i n 

n>0 



Now as k approaches infinity, the sum on the right approaches 
for all convergent (and many divergent) series, leaving, finally, 



1 ^ 

"U a 

n>0 k>0 



2^1 n l2^/ 





(U a - a ), The motivation for this transform should 
■ n n . 

become clear with an example. Ue have 

pi/4 - arctan 1 - 1 - 1/3 + 1/5 - 1/7 + ... 

n 

(-) 

so that a - — — 

n 2n+l 
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n 



n n+1 

1 (-) (-) (-) 

U a ■ - - { + ) » — 

n 2 2n+l 2n+3 (2n+i) (2n+3) 



n n+1 n 

2 1 (-) (-) 2 (-) 

U a - - (— + ) . 

n 2 (2n+l)(2n+3) (2n+3) (2n+5) (2n+l) (2n+3) (2n+S) 



It i3 easy to induce that 

n 
k ■ k! (-) 

U a - , 

n (2n+l) .... (2n+2k+l) 

Setting n to 0, we find 

n 
pi ^-^ (-) l^-^l 2 3 k 

4 Z^2n+1 2^-^3 5 7*" 2k+l 
n > k>8 

The left hand sum has as a ratio of consecutive terms 

' 2n+l 

_ 9 w hich approaches -1 as n grows large. 

2n+3 

To get 9 digits of pi would require about half a billion terms. 
The right hand sum is easier to compute than it looks, since 

k 

its kth term is just times the previous. 

2k+l 

Since this term ratio approaches 1/2 as k grows, only 30 terms 
areneeded for 9 digits (38 bits). This convergence rate of 
1_ bit/ terra is typical of Euler transformed series. 

Note, however, that if by ugly coincidence some 

different series happened to match our pi/4 series for the 

first. 30 terms (and then, for instance, became all zeros), the 

first 30 orphans would have the same numerical 

values, albeit from a different formula, in which 

case few, if any, of the 9 digits would be right. 

Moral j virtually nothing is certain until you have the formula 

for the kth orphan. 
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Generalization of the Euler Transform 

Unfortunately, if the original series converges more rapidly 
than 1 bit/term (term ratio < 1/2), the transformed series 
may converge less rapidly, as we shall see. This is a symptom 
of the fact that Euier's transform is most effective when 
consecutive terms are of nearly equal magnitude and opposite 
sign, so that their averages will be very small, leaving most 
of the sum's value in the orphan (flore orphan than not). ■ 
What we need is an unsymmetr ical , weighted averaging operator 
to achieve this near telescopy in series where each term is 
significantly smaller than its predecessor. We are indebted 
to Rich Schroeppel for suggesting the weights l/(l-r) and 
r/(r-l) be used when the series has a limiting term ratio r. 
This uill actually split the terms of non-alternating series 
(positive r) into pieces with opposite sign. It was an easy 
step from this to realize that the weights need not be the 
same for each term. 

Thus we will denote a pair of weights by s. and 1 - s , 

n n 

and cal I s a splitting- function, 
n 

Ue generalize Euier's transform as follows: 



E-.-E 

n>0 n>0 



(s + 1 - s ) a 

n n n 



- s a 

e a 



*E 

n>8 



(1 - s ) a + s a 

n n n+1 n+1 



which specializes to the first step of Euier's transform when 
s « 1/2. s - leaves the series alone, while s =1 
n ; n ■ n 

translates the series by one term, the zeroth becoming the 
orphan. He wish to choose s , the splitting function, so 

n 

as to make the 

summand both small and simple; small so that upon repeated 
applications of this operation, the value of the sums will 
accrue rapidly in the orphans, and simple so that these 
orphans have a general formula easy to discover and compute. 



PAGE G 



For examples, let s = s, a constant. Transforming 

n 

the geometric 9eries: 



In n n+1 

— - y x ■» s + y d-s) x + s x 

n>0 n>0 



n 
s + (1 - s + s x) N x 



E 



n>0 



2 
s+{l-s + sx)s+(l-s + sx) s + 



E 



k s 

S > (1 - S + S x) = : 

1 - (1 - s + S x) 
k>0 



an identity, except the interval of convergence has moved 
from -1 < x < 1 to 1 - 2/s <x<l or 1<x<1- 2/s, 
depending on the sign of s. 

If we choose s so that 1 - s + s x => 0, all of the terms of this last 

sum vanish except the first, leaving s which is indeed 1/Cl-x). 

Ue note that if s ■ 1/2 (Euler's transform) and 

-1/3 < x < 1/2, "the transformed series converges more 

slowly than the original geometric series. 

Integrating these series from 8 to xs 

n k k 

x ^ — ^ (1 - s + s x) (1 - s) 

- In (1 - x) 



■Er-E 



n>l k>l 

In (s - s x) - In s 
Here, if we again choose s - i/(l-x) 
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. n x k 

£;- - - ElEL 

n>l k>l k 



X 

or - In (1 - x) - i n (i ) 

x - 1 

If x • -1/2, we can compute In (3/2) with an asymptotic 
term ratio of 1/3 using the transformed series, versus -1/2 
for the original. 

If we had chosen s = 2/(2-x) instead of 1/U-x), then 

■1 - s+ s x - - (1 - s ), that is, alternate terms in the 

sum for In (s - s x) - In s will cancel, leaving 

n x 2 k + 1 

x v"^ ( ) 

2 - x 



s-- .-z 



n>l . k>8 2 k + 1 



1 + y x 
or - In (1 - x) » In with y 

1 - y 2 - x 

Here, if x « -1, the term ratio approaches 1/9, yielding In 2 
with almost 1 decimal digit/term instead of 1 bit/term, as 
with Euler's transform. 

The formulas we have seen up to now have been known for a 
long time. Before we get into the more interesting goodies, 
I want to replace the summation notation with one that stress 
term ratios, and reformulate the series transform accordingly. 

R Notation - 



sses 



Let R r ■ stand for the series whose first term is 1 
n>m n 

term is r and whose ratio of r.+lst to nth term is r 

m " . n' 

i.e., 1 + r (1 + r (1 + . . . . 
m m+1 



Thus 



a 
a + a + ...» a R _ n +l . 

m m+1 m n>m a 
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Note that m need not be an integer: 



2n + l 
R 

n>0 2 n +'3 



R 



n 



1 n + 1 

n>- 



Other examples of R notation: 



T- - R - 

jL~J n\ n>l n 



n>B 



In (1 + x) 



£>- 



x 2k + l 
( ) 



k>0 2 k + 1 

2 x 2 k + 1 x 2 
R ( ) 

x + 2 k>0 2 k + 3 x + 2 



F(a+l,b+l;c+l;z) 



c! 
a! b! 



E 

n>0 



(n+a) ! (n+b) ! n 

z 

(n+c) ! n! 



n + a n + b 
n>l n+c n 

(F is the hypergeometric function. Unfortunately, a. b, and c are 
Pictorial ; *™ arenU * to -accomodate the Gamma nota io Tor the 
f actor. al function. The author finds this unnatural and pedantic.) 

uithV^A? 8 * th l S ™ tatl °" ca ™°* express a series starting 
with 0. Also, the (2 k + 1 ) / (2 k + 3) in the logarithm 
formula ,s awkward-suggesting the fol Sowing generalization: 

Let b R r stand for the series 
n n>0 n . . 

■> 

b + r (b + r (b + ... « 

a 0ii 2 

b + r b + r r b +-r r r b + ... 
0B1012B123 
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n+1 

i .e. b R . r 

8 n>8 b n 

n 

b R r ' = 
n n>m n 



Another way to say it is that 



where S satisfies the first order, linear recurrence relation 



S - b + r S 
m mm m+1 

Then our logarithm formula becomes 

2 x 



in (1 + x) 



(2 k + 1) R 

x + 2 ■ k>B 



x + 2 



2 
•) . 



Now our notation is no longer uniquely determined by 
the series — in fact one could cop out and write 



a + a + 
m m+1 



a R 1 . 
n n>m 



The intended convention, however, is that only factors of 
the nth, but not the n+lst term appear to the left of the R. 
(And, of course, the necessary scale factor to make the first 
term come out right.) This will prevent factors of the form 
.(k + a)/(k + a + 1) or its reciprocal from needlessly 
complicating the expression to the right of the R, whose 
primary purpose is to get rid of nth powers and factorials. 

Thus, Ramanujan's amazing pi formula reads 

4 21460 n + 1123 (4 n + 1) (2 n + 1) (4 n + 3) 
— » R 

pi 882 



R 
n>0 



4 3 
128 21 tn + 1) 



(Col lected Papers, US, .eq 39). 



The advantages of this notation a 



re: 
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One can tell at a glance the first term 
and term ratio, which immediately indicates 
the required computation effort and convergence 
rate. The number of digits/term is just 
-log Lim lr I as n grow large, where the base 
n n 

of the log is the radix of the digits. 

The notation is usually more concise, since term 
ratios are usually simpler than terms. This often 
makes it easier to notice when a series is equal 

to, or a special case of, another. 

It suggests a generally more efficient method of 
numerically evaluating n terms of a sum — 
instead of computing the terms and adding them up, 
take r , add b , multiply by r 

n n n-1 

add b , ... multiply by r , and finally add b . 
n-1 

It has suggested two "spigot" methods of 

evaluating sums to unlimited accuracy without multiplying 

or dividing very long numbers. A "spigot method", is 

one that produces digits on demand, computing 

only as much as necessary, yet able to. supply subsequent digits 

upon subsequent demand, without, in effect, starting over. 

It is immediately rewritable as a continued fraction: 

P 



) R 
n n>0 



b + p 
8 8 



+ q 
_0 

b + p 

1- 1 



Although it isn't quite as manipulable as summation 
notation, we can do things- iike . 



dab a -l b 

— x R r x = (k b + a) x R r x 
dx k>8 k k>8 k 

Finally, it simplifies the forthcoming discussion of series 
transformations. 
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Derivation of the Transform 

After applying k-1 splitting functions, thereby extracting k-1 
orphans, we are ready to apply the kth splitting function s 

k, n 

to the remaining f.eries, whose nth term we denote a , so that 

k,n 

a - a , the nth term of the original series. Also, 
0,n n 



k,n 



k , n+1 

a 

k, n 



the nth term ratio of the ktuply accelerated series. Then, 



ra - V~^ a { S + 1 
k > n Z-/ k-.n k,n 



- s ) 
k,n 



n>m 



n>m 



- s a 



+ V^ (1 - s ) a -4- s a 
k,m k,m /^ k,n k,n k,n+l k,n+l 



n>m 



- . s a 
k,m k 



+ . V^ (1 - s + s r . ) a 
» m . Z— 4 k,n k,n+l k,n k,n 



'n>m 



and now if we define u . - 1 - s . + P s 

k *h k,n k,n k,n+l 



Z^ k.n 



s a 

k,m k 



h>m 



+ y u a 
im / / k,n k,n 



n>m 



so s a 



k,m k,m 



is the kth orphan, and the nth term of our new ser 



i es is 



u a • which we define to be a . Thus 
k f n k,n k+ i (n . 



/ -/ k,n 



s a 

k,m k 



n>m 



+ V a 

♦ ni / j k+l,n 



n>m 
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As before, we iterate on k indefinitely, choosing s so - 

k,n 

that the u and thus the a and thus their sums tend to 
k,n k,n 

with k as rapidly as possible. Then nothing will remain of 
the original sum but the aggregate of the orphans. If k is 
ini t i a I l.y j , 



Ea «s a +s a + . . . 
J» n - J * ni j,m j+l,m j+l,m 
n>m 



» s a +s u a +s u u a + 
j»m j,_m j+l;m j,m j.m . j+2,m j+l.ra j.m j,m 

because a - u a by definition. 

k+1 , n ' k , n k , n 

In R notation, this last equation reduces to 

a R r • - a s R u . 

j,m n>m j,n j,m k,m k>j k,m 

Note that we have nowhere used the fact that j is initially 0, or 
even an integer, so that we are free to choose any convenient 
value to denote our original series. " 

Ue now repeat this whole derivation in R notation, skipping fewer 
steps. The reader is advised to understand it, since the main 
results in this paper depend heavily on it.. .. 

For convenience, we assume that the original series 
has been scaled so that a =1, Then 

j.tn 



E a , - R r - (1 - s + s ) R r 
. j» n n > m J.n j,n j,n . n>m j,n 

n>m . . ~ 



- tl - s J R r + s . R r 

j , n n>m j , n j , n n>m j , n 

CI - s ) R r + s + r s R r 

j,n n>m ■ j,n j,m j,m j,n+l n>m j , n+1 



(1 - s ) R r + s + . r s R r 

j,n n>m j, n j,m j,m ),n+l n>m j , n+1 r 



J t n 
(1 - s^ ) R r + s + r s R r 

j , n n>r.i j,n. j , m j,mr j, n+1 n>m j t n 

j,m 

- (1 - s . ) R p +b +r s Rr 

j,n n>ni j,n j t m j,h j,n+l- n>m j,n 

- S + (1 - s + r s ) R r 

j» m j.n jtnj, n+1 n>m j , n 



s + u R r 

J t m j i n n>m j , n 



j,n+l 
3 + u R r 

j » w j , m n>ra u . j , n 

j.n 



s + u R r 

j.m j»m . n>m j+l,n 



u 
k,n+l 

* r ■ r because a - u a ) 

k+1 » n u k,h k+l,n k.n k t n 

k,n 



s + u (s + u R r ) 

i' m j»f" j+lf-m j+l.,m n>m j , n+2 

s + u (s + u . (s . + . 

i' m ■ j.m j+l.fli j+l,m j+2,m 

s R u 

k,m k>j ' k,m 
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r 
j.n 
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Final ly, we have the 



Series Acceleration Formula 



R r 
n>m j , n 



s R u 
k,m k>j k,m 



Our problem then boils down to finding a pair of functions 
s and r which satsify the recurrence 
k,n k f n 



k+l,n 



k,n+l 
r 

u k,n 

k, n 



with u =l_ s +rs 

k.n . K f n k,n k,n+l 

while for some j, r equals the term ratio of the series 
we wish to accelerate. 



Constant Splitting Functions 

Since log and arctan are particularly useful, we will honor 
them with yet another example of a new way to 
get an old resul t; 



- In (1 - x) 



n 
x R X 

n>l n+1 



n + 1/2 2 
arctan x - x R (_ x } 

n>0 n + 3/2 



These are just special cases of the hypergeometr i c function 
with 2nd argument - 1, . 



F(a,l;c;z) 



n + a 
R z . 

n>8 n + z ' 
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Here m - 8 and again, our asymptotic term ratio is z, so if. 
we continue to restrict s to be a constant, l/(l-z) is still 

-k,n 

our best guess to minimize u . This follows from setting 

k, n 

s and s . to s, r to z, and u to B 
k,n k,n+l k,n K,n 

in the definition of u 

k, n 

We choose j « and have 



n + a 



0,n n 



+ c 



1 ■". 1 n + a. 
u « 1 " + 2 

» n 1-z 1-zn + c 



c - a z 
n + c z - 1 



so 



l.n 



n + a 
n + c + 1 



Since this is just r with c+1 replacing c t and c is arbitrary, 

0, n 

we can immediately conclude that 
c+ k - a z 



k,n n + c + k z - 1 



Final ly, 



n + a 
R 2 

n>0 n + c 



k + c - a 



1-z k>0 k + c z - \ 



or more symmetrically, 



n 

- 2 R 2 

n>a n + d 



2-1 n>d n +. a z - 1 



(with d for c-a) , a result usually found by analytic 
methods (See Abramowitz and Stegun, Handbook of 
Mathematical Functions, Formula 15.3.5) . 

(Note that the same transform on this new series yields the 
old one again, so that for every z for which convergence was 
improved, there was the point z/tz-l). for which it was 
degraded, and vice versa.) 

In particular, we have 



arctan x « 



2 

x + 1 



k + 1 
k>8 k "+ 3/2 



2 

x + 1 



the interval of convergence having expanded from + 1 to 



2 2 2 
± infinity. Setting y - x /(x + 1) we get 



arcsin y 

2 1/2 
(1 - y ) 



k 2 

R y 

k>l k + 1/2 



k>8 ( 



2 k + 1 

i 

+ 1/2 
) 



where 



J j! 
( ) » 

k k! (j-k)! 



which, for integer k 

" j (j-D.- 



J 

( ) 

j-k 



(j-k+1) 



k! 
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so that 



n - 1/2 

( ) 

k 



2 n 2 k 
( ) ( ) 

2 k k 

n k 
( ) 4 
k 



Integrating the arcsin equation: 



arcs in y 
2 



2 k + 2 



£ _H 

k + 1/2 
k>0 (2 k + 2) ( ) 

k 
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2 k>l • (k + 1/2) (k + 1) 



y . k 2 

:-- R u 

2k k>l k + 1/2 



This result, although known, is peculiar in the light of 



arcsin y 



k - 1/2 2 k + 1 

E ..J. . 

k>e 2 k + i 



y k + 1/2 2 

R y 

2 k + 1 k>8 k + 1 



wherein the coefficients are nearly the reciprocals 
of those m the previous sum. The relationship 
is even more striking if we write 



k 2 

y 



arcsin y - R 

2k 1 k + 1/2 
k>- 

2 



which says that. if we offset the index of summation by 



-1/2 in the formula for arcsin y, u 



e get y arcsin y. 
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Other Sp! i tt ing' Functions 



Actually, for all but the simplest term ratios, most splitting 

U !^'. C ° nSt ? nt or.otheruise. wi I I lead to prohibitively 
complicated u 's after just a fewk. 
k, n 

So the main problem is to find an s such that u 

k.n k,m 

stays simple. The easiest way 

to do this seems to be to start with r and choose the s 

0,n e,n 

which results in the simplest r you can find. Then repeat 

l,n 

the procedure for r , seeking the 9 leading to the simplest 

l.n . - l,n 

possible r .With some skill and luck, you find a ' 
2,n 

sequence s , s , ... which produces a sequence of u's 
o f n l,n 

whose general formula is di scernable. . 
For example, let us attack 



zeta 2 



1 1 1 

— + — + — + ...■- R 
2 2 2 

12 3. 



n>l 2 
(n+1) 



Formula "anri^Tl" "'""* *° Ch °° Se j ' 1 ln the Acceleration 
rormuia, and so define 



l,n 2 

(n + 1) 



Ue certainly can' t try s 



y s - Lim l/(l-r ), since the asymptotic 
l»n n n 



term ratio is 1, so let us try a ver^ simple function of n 



s ■ a n + b 

l,n ■ 



with a and b to be determined so as to simplify 
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2 

U-a) n + (2-a-2b.) n + 1-b 

■ u . - 

l,n 2 

(n + 1> 

A good heuristic is to reduce the degree of u, which we 
can do twice by choosing, a - 1. b - 1/2, so that 

1 
u = 

l,n . 2 

2 (n + 1) 

Then 

2 
n 

r - , 

2,n 2 

(n +-2J 

Trying the same form again, 

s ■ - a n + b, 
2,n 

we have 

2 

(l-3a) n + (4-4a-4b) n + 4-4b 
u - 

2,n 2 

(n + 2) 

Analogously, we choose a « 1/3, b - 2/3, so that 

4 

U n _■ 

2,n 2 * 

3 (n + 2) 

Then 

2 . . 
n 

r m 

3,n 2 

(n + 3) 

which suggests 

2 

n 
r » 

k»n ' 2 

(n + k) 
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St ick ing u i th 



s -- a n + b, 

k, n . 



2 2 2 

tl-a(2k-D) n + (2k-ak -2bk) n + (l-b)k 

u ■ '■ 

k.n 2 

(n + k) 



which says 

1 3 k - 2 

a . , b » 

2 k - 1 2 (2 k - 1 

3 
k 

u - . 

k,n 2 

2 (2 k - 1) (n + k) 

and, sure enough, 

2 
n 

r - , 

k+l,n 2 

(n + k + 1} 

So final I y, 

2 n + 3 k - 2 



s - 

k,n 2 (2 k - 1) 



and 



3 
3k; k 
zeta 2 = R 

2 (2 k - 1) k>l 2 

2 (2k- 1) (k + 1) 



3 . k 1 

R . 

2 k k>l (k + 1/2) 4 



with a convergence rate of 2 bits/term. 
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D - • 2 

But this is just a special case of our arcsin formula 
wi th y - 1/2, so that 

2 
3 8 2 1 pi 
zeta 2 - - - arcsin - » . 

'.22 2 G 

(I am indebted to Dr. John U. Urench, Jr. of the Naval Ship 
Research and Development Center for pointing out that this 
zeta 2 formula is derived as a Markov transformation in 
Section 156 of Knopp' s Theory and Application of Infinite 
Series. Later, however, we will improve this to a B bits/term 
formula, which has a better chance of being undiscovered.) 

This/example has nicely illustrated the actual practice of the 
method: 

■Acceleration Recipe 

1). Urite the term ratio. r 

0,n '.- 

2) Guess a form for a splitting function s , cunningly 

B,n 

including some undetermined coefficients. It probably should 
be asymptotic to l/(l-r ). 

B,n 

3) Choose these coefficients to simplifly 

u «=l_ s +s r 

0." 0,n 0,n+l 0,n 

u 

0,n+l 

4) Wri te r - r 

l.n u. B, n 

0,n 

5) Guess r so that it agrees with r and r 

k ' n 0,n l,n 



u ' 
k.n+1 

5) Writer * r and determine the 

k+l*n u k,n 

k,n 

coefficients in s to achieve agreement with r 

k ' n k,n 
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This ail boils down to searching for rational solutions to 
the recurrence relations defining u and r 

k, n k,n . 

and there may be much to gain from exploring the literature of 
nonlinear difference equations, which I have as yet neglected. 

Now let's try harmonic numbers. 

Ill 1 p 
h . » _+_+_+..,+_ „ R 

P 1 2 3 p n>l n+1 



We. have not equipped our transform for finite series — it 
would evidently produce a sum of two series of orphans, one 
from the left and one from the right. Ue can, however, use 
the transform to see what the one on the left will be. 

I f we agai n try 

s - a n + b 
l,n 

with 



r - 

l,n n+1 



we f i nd 



n - b + 1 

U a 

l,n n+1 



and 



n 

r a 

2,n n+1 



independent of a and b! This means that we can extract 
an arbitrary series of orphans from the left without 
affecting the term ratio at all.. 
Other-choices of splitting functions are equally 
unenl ightening. So it wants an infinite series? Ue 
g i ve i t an infinite series. 
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111. 
h - _+_+_ +#># 

P 12 3 

11 1 
- ( + + + tm t 

p+1 p+2 p+3 



E--- L -E— 

^— ' n n+p *— ' n (n + 
n>l n>l 



A slightly less dubious way to get this is to let m approach 
infinity in the cute identity 

m p 

\""* P ^-^ m 

^-— ' n (n + p) ^^n (n + m) 

n»l n «l ' 

but this seems to rely on p being an integer— a restriction 
which we will profitably waive in-order to conquer the di gamma 
(Pst) function while we are at this: 



h " 7 » Psi (1 + p) + enigma 

P- ^— ' n (n + p) 

n>l 

where enigma is Euler's constant - .5772156649. . . and is usually 
spelled gamma. Ue have 

n tn + p) 

r -. — . : 

l,n (n + 1) ( n + p + i). 
Using good old . 

s - a n + b", . 
l.n 



we get 
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l.n 



Cl-a) n + (p-2b-a+2) n + (1-bMp+l) 
(n + 1) (n + p + 1) 



which suggests 



a « 1 and b 



P + 1 



so that 



and 



(p - 1) (p + 1) 

j . 

l.n 2 (n + 1) (n + p+. 1) 



n (n ■+ p) 

r * — 

2 f n (n + 2) (n + p + 2) 

Let us leap to the conclusion that 

• n (n ■+ p) 



k,n (n + k) (n + p + k) 



Then 



u 



k,n 
so that if 



2 2 

(l-(2k-l)a) n + (p-2kb-(k +jk-p)a+2k) n + (l-b)k(p+k) 

• (n + 1) (n + p + 1) 



giving 



1 3 k + p - 2 ■ ■ 

a » and b » 

2 k - 1 4 k - 2 



2n + 3k + p-2 



k,n 2 (2 k - 1) 



then 



k (k + p) (k - p) 

J m ' 

k,n 2 (2 k - 1) (n + k) (n + k + p) 
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and lo and behold 

n (n + p) 
r * 

k+l,n (n + k + 1) (n + p + k + 1) 



Final I y 

h 

. P p + 1 n>l n + 1 n + p + 1 



P n n + p 



P ■ .P 
R r - s R u 

p+1 n>l l,n p +1 k.l k>l . k,l 



P 3 k + p . k (k + p) {k - p) 
, p ; 

p + 1 2 (2. k - 1) k>l 2 (2 k - 1) (k + 1) (k + p 4- 1) 



P (3 k + p) k - p 1 

R ___ 

2k (k + p) k>l k + 1/2 4 



k - p - 1 
(3 k + p) ( ) 

.-£-._ !!_'. 

*-* 2 k 

k>l k (k + p) { ) 

■ k 

As a consistency check, we can divide through by p,- then 

let p approach- 8, yielding the zeta 2 equation we got earlier 

Note that if p is an integer > B, the k-p in the numerator 
of the ratio wi I I . terminate this series on the jth term, 
Ue have our finite series back, only i t converges faster! 

For: instance, 

111111. 1 761 
1 +- + - + - + . + _ + - - + __. m 

2 3 4 5 G 7 8 2S0 



44 49 238 5-46 13 23 1 
3 15 165 12 585 1386 45845 51480 
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Now if we are computing harmonic numbers as exact nationals 
it is easier to add up the first series than the second. 

But if we want big ones, this becomes impractical. 

(h^ involves 83 digits.) If we resort to decimal approximations, 

"• -can still sum the harmonic series using approximate reciprocals, 

but even this grows impractical as p gets to be a few hundred. 

There is a divergent series involving log p, Euler's constant, 

and Bernoulli numbers, but even with an abundance of digits of 

these, you cannot extract many more than 2 pi p/ln 10 digits 

The formula we have just derived seems 

to solve our problem to the tune of 2 bits/term, but on closer 

inspection, we see that the power of p in the numerator 

of the term ratio is one greater than in the denominator, 

so that if p is large, the series will not start to converge 

until k reaches (p - 2)/5. 

To get the most out of this formula, we should notice a degree ■ 
of freedom in the series transform that we have not yet 
exploited— namely, that setting n to m in the s and u of 

k,n k,n 

the transformed series corresponds to starting the original 

series at the mth term. (He must also remember to scale both 

series so that the mth term of the original series is of the 
desired size.) Thus we have 

11 1 : ■ 
"lip m ■ m+1 m+2 

111- 
- { — + _„_. + + ... ' 

m+p m+p+1 m+p+2 

- h - h 

m+p-1 m-1 



1 1 



■ E - - - : - - T ----- 

^—f n n+p Z— 4 n (n + 
n >rc n>m 



P) 
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p (3 k + 2 m + p - 2) k (k + p) (k - p) 
R 

2 m (m + p) (2 k- 1) k>l 2 (2 k - 1) (k + m) (k + m + p) 

p (3 k + 2 m + p - 2) k (k + p) (k - p) 1 
R _ ■ 

2 m {m + P } k>i (k + 1/2) (k + m) (k + m + p) 4 



■ E 



k + p k - p 
p(3k + 2m + p + l)( }( ) 

k k 

_ 2kk + mk + m + p 

k>8 2 m (m + p) (2 k + 1) ( ) ( ) ( j 

. k k k 

Now if m is about as large as p, we see that this new term ratio 
is about 1/8, so that we have an almost 1 digit/term method 
for summing the reciprocals between half a trill ion' and a trillion, 
for instance. But' we can. use this same trick for those between a 
quarter tri I I ion and hal f a tr i I I ion, etc., so we really 
can accurately compute very large harmonic numbers subject to 
a cost factor of log p. 

2 

This example shows that once we successfully accelerate a series, 
we have a handy way to sum finite but large intervals of that 
series, since the starting value of the summation index is merely 
a parameter in the new series. In fact, it need not be an integer: 

Ue know 

1 ^ 2z ' 

cot z - -- +y^ — . — _„. 

z /-/ 2 2-2 

k>l z - k pi 

then 

pi 2 COt pi Z • - 2 h - 1 - 

. . -z,2z 



k + 2 z k ~ 2 z 
(3 k + 1) ( ) ( .) 

-.i + Y ... ... -._ t..._ 

*—* k + z k - z 2 k 

k>3 (2 k + 1) ( )(■'.)( ) 

k k . k 



2 2 
k (k - 4 z ) 
-1 + (3 k - 2) R . 

k>l 2 2 

(4 k + 2) (k - z ) 
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a 2 bits/term cotangent formula. 
Similarly for Catalan's constant: 

1111 

2 2 2 2 
1 3 5 7 . 



■ 1 1 1 

' 2 (— + -- + _- + 

2 2 2 

1 5 9 

1 1 1 

-(-- + „ + „ + 

2 2 2 

1 3 5 

h h 

1/4, p 1/2, p 

p->0 IS p 4 p 



2 

l v x pi 

".2 ^ 2 " T 
k>B k + 1/4 2k 
( ) ( ) 
k k 







3 


2 


3 




(k + 1) 


pi 


- 


R 


. 




2 


k>B 


2. 
(k + B/4) (4 k + 2) 


8 



Later on, when we accelerate eta 2, we will get a single 
series formula for Catalan's constant. 

The next example will emphasize a very important point. Ue 
recall that our original goal was to accelerate the particular 
series - 

1 . 1. "' 1. 
^_ + + + 

1 (1 + p) 2 (2 + p) 3 (3 + p) 
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with 

n (n + p) 
r - . 

l.n (n + 1) { n + p + 1) 

but in the process, we accelerated the more general form 

n (n + p) 
r - 

k,n (n + k) (n . + p .+ k) 

This means that if we wish to accelerate some other series, and 
its term ratio is expressible as 

r 
a,n+b 

for some a and b, then we can use the transform we already found, 
summing from k - a instead of 1, and setting m (the initial value 
of n in the original series) to b+1 instead of 1, Thus. if we 
wish to accelerate the apparently very different series 

4 2 -12 -112 -113 2 -113 5 2 
— - 1 + („) + ( } + { ) + ( ) 

pi 2 2 4. 2 4 S 24B8 

(Jol ley's Summation of Series, formula #274) 

2 n - 1 2 
we find the term ratio » ( ) 

2 n + 2 

But this is just r with p = 8: 

3/2,n-l/2 

.2 
2 n - ■ 1 2 
R ( ) „ R _______ 

n>0 2 n + 2 n>-l/2 2 

(n + 3/2) 

3 
■3 k - 3 k 
_ R 

2 (2 k - 1) 3 3 

k>- (2 k - 1) 
2 

■3 (2 k - 1) ( 2 k +"1) ' 
- _ R 

8 k>l- 2 

32 k ik- +1) 
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so 

2 2 k 3 
(2 k - 1) k ( } 

P ! ^ 8 k 
*>1 2 




t^ ^ho'fir^ f 6 ° r ; ^'f^ PaVeMe ° f P -«ption Technology 
Inc., who first suggested transforming the harmonic series. 

Let us now return to 

pi/4 - arctan 1 - 1 - 1/3 + 1/5 - 1/7 + ... 
to see how we can improve on Euler's transform. Ue have 

2 n + 1 

r m- 

0,n 2 n + 3 

l/2 Ch A^PrT" 63 7 1 ' SU ^ estin 9 a splitting function approximating 
1/2. A very s.mple non constant such function is apparently 

n + a 
s - 

8,n 2 n + b 

with a. -and b chosen to simplify 

2 2 2 

------!!!i_]„!j 2b +10b - 4ab -12a+4.) n + 3b + 5b-4ab~6a 

* n (2 n + 3) (2 n + b7~"(2~n"~b~+"27"~~~~~~ 

Reducing the degree of the numerator as before, we let 
D = ^ a - l. Then 



u 



2 
2 <(2 a - 2)-'n + 2 a - 2 a - 1 



. B ' n (2n + 3) (2 n + 2 a - 1) I2~nl~2~l~71) 

Instead of simply taking the obvious choice of a - 1 
let us proceed directly to 
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r 



u • 
B,n+1 

1 1 n u ■ B, n 

8,n 

n 

^ 2^n + 1 2 n + 2 a - 1 (2 a - 2) n + 2 ^ - 3 

2n + 5 2n + 2a + 3 2 

(2 a - 2) n + 2 a - 2 a - 1 

to see which choices of a allow at feast one of the denominator 
terms to cancel a numerator term. For each pairing of 
a term from the denominator with one from the numerator, we 
scale them so that the coefficients of n are equal, then solve 

maPcvma a *u Ch maKeS the Wh0le terms ec * ua! - Uhen us5n 9 
HALbYHA, the very powerful symbolic math program at HIT 

?u S b?Sr,. P ?AK,? rm thiS °P eration succinctly by invoking 

the KhSULTANT function,, whose existence and usefulness 

were revealed to me by Professor Hal Abel son. The resultant of 

the numerator and denominator of the above expression is 

IB 2 2 2 
" 2 a (a " 1} , {a ~ 2) (a - 3) (a + 1) (2 a - 1) (2 a - 3) 

and in general, the resultant of two polynomials 
is a polynomial in their coefficients which is 8 when they have a 
common factor or" their leading coefficients simultaneously 
vanish. They are painful to compute, but will no doubt come 
oack into popularity now that machi nes wi I I do them for us. 

The series in question is so easy to accelerate that any of these 
choices of a will lead to a ratio which is in turn easy to split, 
etc., and in this way we can generate dozens of pi formulas of 
varying convergence rate. Perhaps the most interesting, yet easily 
found formula lies along the path indicated by choosing a - 3/2 
I nen s 

(2 n + 2) (2 n'+ 3) 
r- w — : ■ 

!» n (2 n + 5) (2 n + 6) 

The simplest guess of r consistent 

k,n 

with k - and 1 is apparently 

(2 n + k + 1) (2 n + k + 2) 

r . - ■ 

K ' n . (2 n + -3 k +2) (2 n + 3 k + 3) . 
which is, in fact, supported by 
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2 n + 5 k + 3 

s - ' 

k,n 4 n + S k + 4 

(which can be found with some pain by determining a and b 

n + a : ■ 

m s = ~ '" so that r conforms to our conjectured 
k t n 2 n + b k+l,n 

formula for r ). Then 
k,n 

(k + 1) (2 k ■+ 1) (2 n + k + 1) 

u •« 

k,n (2 n + 3 k + 2) (.2 n + 3 k + 3) (2 n + 3 k + 4) 



so that 



P* 1 5 k + 3 (k + 1/2) (k + 1) 2 
— , R — „ 

4 4 k>0 (k + 4/3) (k +5/3) 27 



E 



(5 k + 3) (2 k) ! k! 



k+1 
*>£ 2 (3 k + 2) ! 






at better than 1 decimal digit/term. Interestingly, thi 
—rtes.can be simplified by rerunning the transform just 
e step on this new series, with 



se 

one 



(5 k + 8) (k + 1/2) (k + 1) 2 

* _ 

0,k (5k+ 3) (k + 4/3) (k + 5/3) 27 



and 



6 k + 4 



0,k 5 k + 3 

which happens to be the reciprocal of the old s 
Then 



0,k 



so 



B (k + 1) (k + 2) 
3 (3 k + 4) (5 k + 3) 
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(k + 1/2) (k + 3) 2 

r m ■ 1 

1,K (k + 5/3) (k + 7/3) 27 

and thus 

P' 1 '5 (k + 1/2). (k + 3) 2 

— - 1 R 

4 24 k>8 . (k + 5/3) (k + 7/3) 27 

or 

. 

k 3 k + 4 

k>0 2 (2 k + 1) ( ) 

•k + 2 

But this formula is just what we would have gotten had we 
accelerated the series 

P' 1111 
■ -"'■- 1 - — + + - . . . 

4 3 5 7 3 

and then restored the Initial 1; which is to say that 

s R n u can be simplified further than s. R u . 
k ' X k ^ ka k.B k>0 k,0 

This suggests that when we get an acceleration formula, we 
should check to see if -there are any preferable starting 
values of the summation index which differ by an integer 
from the the one initially intended. The brute force way to do 
tm s is to rewr i te ■ 

s 
■ _ k,m+l 

s R u as s R ^__ u 

k,m k>0 k,m . .' 8,m k>0 s k,m 

k,m 

^ T*" - ak ? the / e f ultant -with respect to k of the numerator 
and denominator of the expression, on the right to determine 
the values of m for which it reduces. 

Actual ly," the fact that there were two ways to find that 

bT a npr ^JV"! f^egoing: pi/4, formula was no coincidence," 
but rather a manifestation of a beautiful duality between a 
series and its transform: a 
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If 9 ■ splits the series whose ratio is r to make the series 
k v n k,n 

s 

■k+l.n 
with ratio r u , then s - 1/s splits 

n » k s k.n. n,k k,n 
k,n 

r back to r . (Note the switching subscripts.) 
n,k k,n 

This can be verified directly from the recurrences defining 
r and u 
k f n k,n 

Thus that extra split by the reciprocal of s that we performed 

k,0 

on the transformed series was equivalent to the orphaning of 
the first term of the original series via the dual; transform. 

Just to make sure the reader understands this duality, an example: 
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As we will iater show, 

n + 3k + l 

s - 

k f n 2 n + 4 k + 1 

sp 1 its 

2 n - 1 
r » 

k,n 2 n + 4 k + 1 

to make 

2 k + 1 2 k + 2 

u ■ 

k,n 2n + 4k + l 2 n + 4 -k + 3 

so that the transformed series will have ratio 

3 k + n + 4 . 2 k + 1 2 k + 2 

r m . ._ 

n,k 3k + n + l 4 k + 2 n + 3 4 k + 2 n + 5 ' " 

By dual i ty then, 

4 n + 2 k + 1 
s - — 

k t n 3 n + k + 1 

spl i ts 

3n + k + 42n + l 2n + 2 

r - , ; 

■k-.n. 3 n + k + 1 4n + 2k + 34n + 2k + 5 

to make 

■ 3 n + k + 2 2k-l 

U m 

k,n 3n + k + 14n + 2k + 3 

with the ratio 

2k-l 
r - : — - ' # 

n.k 2 k + 4 n + 1 



Zeta 3 
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111 1 

zeta 3 * — + — + — + ... - K~^ -- 

3 3 3 2-j 3 

12 3 n >l n 



s 



k-1 



2. <<— ' 3 2k 
k>l k ( ) 

k 



the improvement being asymptotical I y 2 bits/term versus 8 
bits/term. To be honest, I first got this formula from a 
Kummer. transform (omitted for brevity), so it is probably known. 
Ue have. 



B.n 3 

(n + 1) 



I f we choose s 



0,n 



a n + b, we find that a must be 1/2 t< 



make the degree of the numerator of u less than the degree of 

0,n 

the denominator. But then we will find that there is, no value of 
b for which r is again 

l,n 

a ratio of cubics, and experience indicates that the larger ■ 
the degree of the term ratio, the harder it is to find a splitting 
function which keeps the degree from increasing yet further. 
The next simplest splitting function asymptotic to n/2 is a 
quadratic divided by a linear polynomial, which will 
contain 3 coefficients to be determined. Now we could begin a 
massive search through al ! values of these coefficients for which 
r simplifies, but first, let me suggest two heuristics useful 

1 ♦ n 



for limiting such searches — 
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1) Choose the denominator of the splitting function s to have 
factors in common with the denominator of the term ratio. 

2) Choose the numerator of 1 - s to have factors in common with 
the numerator of the term ratio, 

Uhile these heuristics are not foolproof, a glance at the formula 
defining u will show why they tend to produce simple u's. 

k,n 

If ue use these heuristics to eliminate all three bugger 
factors, we get 

2 
n 
s - 1 + 

B,n 2 Cn + 1} 
which rewards us with 

2 
■ n (n + 1) 

r - « 

■ ■ l.n 2 

(n + 2) (n + 3) 

suggesting 

2 

n (n + k) 

. r * — ' \ 

. k.n 2 * 

(n + k + 1) (n + 2 k + 1) 

But there is no reason to believe that s should be asymptotic 

k,n 

to n/2 for all k. To find out what it should be, we 
temporarily set s to an and find that a should be 
k, n 

1/(2 k + 2) in order to reduce the degree of u . Now the 

k,n 

factors of the numerator and denominator of r are not all the 
same, so which do we choose for s? Ue need another heuristic; 

3) Choose the largest factors. 

This dictates 

. . - 2 
(n + k) 
s = i +. 

k.n .2 lk+1) (n+2k+l) 
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which is indeed the winner among the eight candidates. Then 

2 2 

(k + 1) (n + k) 

u ■ « . ' . 

k,n 2 

tn +. 2 k + 2) (n + k + 1) (n + 2 k + 1) 

and since s « 5/4, 
k',1 



5 - k 

zeta 3 * R 

2 k>l 2 (2 k + 1) 

4 k 



Which is the R way of saying 



k-1 - 

. ! E ...... . 

2 ^— ' 3 2k 

k>l k ( ) 

k 

• 2 

Again we notice a similarity to the arcsin formula, this time 

with y « i/2, but with with an extra factor of 1/k in each 

term, which we can explain by dividing by y and integrating: 



2 
f\ (asinh y) 
. zeta 3 » 10 / dy 

Jo y 

where asinh is inverse hyperbolic sine. 

In an effort to rewrite this as the usual integral' definition of 

zeta 3, I tried substituting 

y - ' sinh t/2 . . 

2 - 
Use coth t/2 - 1 + 

t 
e - 1 
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and 



f(x) d * 'J Q fM dx - / f( x+a j dx 

and 

n-1 

.1 x t 
zeta(n;x) - ; ^ 

(n - 1) ! t 



e - x 



where 

k • 

zeta(njx) - \ — f 

k>i k 



the generating function of the -nth powers. (zeta(n;l) =. zeta n 
and zeta(l;x) - - log (1 - x)}. Then expanding the numerator 
into 3 terms and using. 

2 
1 pi 2 

zeta(2; — ) {| n $) 

2 15 

* • 

(see J. Edwards' Treatise on Integral Calculus, p, 285+) , 
one winds up wi th 

2 
5 1 pi 5 3 
zeta 3 - - zeta (3; — ) + — log <J>- - - ( log *) 
4 2 6 B 
*■ 

where 6 is the golden ratio 

2 asinh 1/2 
- 1.618... - * - 1 = e 



i i i 



Surprisingly, this formula appears in problem XXVI I 41 ,,, 

of Edwards, credited to the ingenious manipulations of Landen. 

Landen and Ramanujan independently found the analogous relation: 

2 
8 . 1 2 pi i 3 

zeta 3 = - zeta(3; -) + i og 2 (log 2) , 

7 2 84 42 
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which suggests that there is at least one more of these identities, 

2 
with ib instead of 4 or- 2. Perhaps this "is what was garbled to 
produce the obviously false equation in problem 41 iv. 

Clausen, in Cre lie's Journal, vol. 5, (1838), p38, gives 

3- n 



pi n 

zeta 3 - — — R ; — 

2. 1 2 (2 n + 1 
8 n n>- 
2 



R 



8 n 



2 n>l 2 (2 n + 1) 



Noting that these R's are integrals of (arcsin y)/y and 



(arcsin y) /y respectively, we can perform the analogous 
manipulations (starting with y = sin t/2i) to get 



[->..■■ 5 pi 13 zeta 3 

3 3/2 36 
n>8 (3 n + 11 2 3 



E 



and 



E 



(-) 



5 pi 13 zeta 3 

9/2 3B 



n>8 (3 n + 2) ' 2 3 
by separate consideration of the real and imaginary parts. 
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Returning to the zeta 3 acceleration formula, we note that setting 
n to 1/4 instead of 1 gives us a formula for 



1 1 1-1 

— + — + __ + „_ + ttt 

3:3 3 3 

1.5 9 13 



11 1 1 . 

- { - ■+ - + _. + ... 

2 3 3 3 

1 3 5 

1 11 

+ — - — + ---.,. 

3 3 3 

1 3 5 



1 zeta. 3 pi 

- (zeta 3 ~ + ' 

2 8 32 



so 



- IS (k +. 1) 



3 2 
7 zeta 3 pi 88 k + 112 k ■+ 41 
— + R 

15 S4 2 k>e (8 k + 9) (8 k + 13) 

48 (k + 1) (4 k + 1) 
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Splitting Function Addition Formula 

Sometimes it is desirable or even necessary to choose a sequence 
of splitting functions whose formulas alternate between two different 
forms. The most useful special case of this is when 

every others is identically 1. This has the effect of orphan.} ng the 
first term of each partially accelerated series, replacing n by n + 1 
in -the term ratio. As we shall see, this is useful for increasing 
convergence rates, but more important, it allows us to use 
acceleration formulas which would otherwise involve division by 8. 

Consider again the problem of accelerating zeta 3, with 

3 

n 

r - , 

B.n 3 

(n + 1) 

Before finding better heuristics, I found, by analogy with the 
Kummer transform, 

2 ' ' 

n 

S . ■ — - r 

0,n 2 (n - 1) 

giving 

■ 1 

u - 

0,n (n-1) (n+1) 

.2 

n (n - 1) 

and r - 

l,n 2 

(n + 1) (n + 2) 

But s and u don't exist, a problem which we might sidestep bu 
0,1 . 0,1 

remembering to sum from n « 2 instead of 1, adding back in the 
initial term when we finish. .Blithely proceeding, we wi I 1 eventually 
di scover the pair 

2. 

n (n - k) 

r - ' . 

K,n 2 

(n + 1). (n + k + 1) 

2- 

n - - 

k,n - 2 (k + 1) (n - k - 1) 
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with 

2 
(k + 1) 

u - - : t ' 

k,n (n + k + 1) £n.- k - 1) 

But this says that if we start with n-2, we will get a zero 
denominator in the k-1 orphan. So we must translate the 
k-1 series by a term as we did for k - 0, In fact, ue will 
have to do this for every k. A nice way to think of this 
is that before the application of the 

kth splitting function s , we interject the splitting function 

k,n 

s - 1. Thus, before the k*8 split, r. will become 
n B,n 

3 
. (n + 1) 



3 
(n + 2) 



necessitating the replacement of n by n + 1 in s , which 

B.n 

evades the division by 8. By extension,, r will become 

k,n 

.2 
n (n + k) 



2 

(n + k + 1) (n +.2 k + 1) 



just before the k + 1st replacement of n by n + 1, 
after which the value of s must be 

■ ■ k , n 

2 
■ (n + k + 1) 



2 (k + 1) n 



(i.e. n + k + 1 replaces n) . . Ue note that the r formula now 
corresponds to the one we got using the heuristics, but the s 
formula does not. This, is because the s in the earlier acceleration 
combined the effects of the current s and the s - 1 preceeding it. 
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How could we have combined these two alternating spitting function! 
to get the earlier one, had we not already known it? "Luckily", 
some straightforward algebra solves the more general problem of' 
finding the sp I i tting- function s equivalent to successively 

n 

splitting by p and then q : 
n n 



Splitting Function Addi t ion formula 



P + q (1 - "p + p r) 

n n n n+1 n 



In our case, p « 1 and 

ri . ' 

.2 
(n + k + 1) 

q » s » , 

n k,n 2 (k + 1) n .. 

Using 2 

(n + k) n 

n k,n+k 2 

. tn + k + 1) (n + 2 k + 1) 

our new, compound splitting function becomes 

2 
(n + k) 
s - 1 + 

k,n 2 tk + 1) (n + 2 k + 1) 

in accordance with the earlier transform. ■' ■ 

The main point of all this is that once ue successfully transform 

a series by finding an s-r pair, we can find any number of new s-r pairs 

by interspersing the unit splitting function using this important 

special case of the addition formula. Ue summarize it under the 

name 
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Translation Transfon 



s <- 1 -. + s r 

k » n k,n+k+l k,n+k 



r <-■ r 

k,n k,n+k 



Rich Schroeppel first proposed the technique of periodically 
translating a series while accelerating it, originally for the 
purpose of improving the convergence rate. Let us use this on 
some accelerations to see how this improves on the improvement, 
First we derive another pi formula from arctan 1: 



"0,n 



2 n - 1 
2 n +.1 



(uith n starting at 1). Heur istical ly choosing the 
same denominator for the splitting function, while 
requiring it to approach 1/2 for large n, we get the form 



0,n 



n + a 



2n + l 



giving 

2 (C2 - 2 a) n + 2 - a) 
B,n (2 n + 1) (2 n + 3) 

suggesting a - 1 so that 

2 n - 1 

r « .' 

l.n 2 n + 5 

suggesting 

2- n - 1 



k,n 



2 n + 4 k + 1 



requiring a « 3 k + 1 (in the formula for u not shown) 

k , n 
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so that 



g i v i ng 



n + 3 K + 1 
s . . , — 

k,n 2 n + 4 k + 1 



(2k+ 1) (2k + 2) 

u - : 

k.n (2 n + 4 k + 1) . (2 n + 4 k + 3) 



Then 



Pi 3 k + 2 (k + 1/2) (k + 1) 1 

; R . 

4 3 k>0 (k + 5/4) (k + 7/4) 4 

with a modest 2 bits/term. (The author is mildly embarrassed to later 
discover that this formula is merely the result of pairuise grouping 
the terms of Euler's transform, of the arctan 1 series, and could have 
been computed directly by using p « 1/2 * q .in the addition formula.) 

n . h 

Using the Translation Transform, ue get 

2 n + 2 k - 1 

r » ; 

■k,n 2-n + Sk + l- 

n + 4k + 2 2n + 2k-l 
. s - 1 - 

k »n - 2n + Gk + 12n + Gk + 3 

These gi ve 

pi 2 (7 k +-B) k + 1/2 k + 1 1 

.4 15 k>8 . - k + 7/G k + 11/6 27 

with a healthy log 27 - 1,43 decimal digits/term. Ue could have 
gotten this pi formula directly by accelerating arcs in 1/2, for which 



r 



(2 n + 1) (2 n + 2 k + 1) 

k,n 4 (2n + 2k + 2) (2 n + 6 k + 3) 

(2 n + 2 k + 1) (2 n + 2 k + 3) 
s « 1 + 

*>" "' 3 (2 n + S k + 3) (2 n + £ k + 51 



and 
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2 
2 (2 k + 1) (k + 1) (2 n + 2 k + 1) 

u - - ; , ■ 

k » n ' (n + k + 1) (2 n + G k + 3) (2 n.+ G k + 5) (2 n + G k + 7) 
starting n and k at 8. 

Recalling the fact that offsetting n by 1/2 changes arcsin y 

2 
into (arcsin y)/y, we have 

2 2 

pi 7 k + 5 - (k + 1) 

18 .3 (2k+l) k>0 3 (3 k + 4) (3 k + 5) 

k . 3 
2 (-) (7 k + 5) k! 

9 JLd (2 k + 1) (3 k + 2)! 

k>e 

This last formula can also be found by a straightforward acceleration 
and translation transform of eta 2: 



2 

1-111 pi 

eta 2 - — --- + — - — .... __ 

2 2. 2 2 12 ■" 

12 34- 



n ' (n + k) . 
Here r » 

^»n (n .+ k + ■ 1) (n + 2 k + 1} 



n + k 
and s - 1 

K»n 2 (n + 2 k + 1) 
(consistent with the heurisitics) 

2 

(k + 1) (n + k)' 

so u « '. 

k,n (n; + k + 1) (n + 2 k + 1) In + 2 k + 2! 



2 

3 k 

and eta 2 » - R ■— -^ ; 

4 k>l (2'k + 1) (2 k + 2) 
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(one half times Mdrkov' s zeta 2 formula). 
Then, translation transforming; 

n + k n + 2 k 

r * . 

k.n n + 2 k + 1 n + 3 k + 1 



(n + k) (n + 2 k) (n + 4 k + 3) 
s - 1 



k » n 2 (n + 2 k + 1) (n + 3 k + 1) (n + 3 k + 2) 



2 
(k + 1) (n + k) (n + 2 k) 



u 



k ; n (n + 2 k + 2) (n + 3 k + 1) (n + 3 k + 2) (n + 3 k + 3) 

whereupon s R u gives us 3/2 times the earlier formula 
0,1 k>0 k,l 

2 

for pi /18. Curiously, s factors, so that 

k,5 

, 2 k 3 

P' 1. 1 1 1 (-) (k + 3) (7 k + 12) (k + D! 

12 1 4 9 + lG " " 2j 2 " " * 

k>l k (2 k + 3) (3 k + 4) ! 

Repeating the translation transform on the pi/4 formula, 
we get a series with 6 bits/term, but these terms have an 
unfavorable quadratic in their numerators. Noting, however," 
* hat r - 's the term ratio for In 2, we'qef 
. . 0,3/2 

In 2 - s R u 

k,3/2 k>8 k,3/2 

14 k + 11 (k +1) (2k +1) 
- R — _ 

IB k>0 8 (4 k + 5) (4 k + 7) 



3 (14 k + 1.1) 



/:: --- 4 k + 1 k+2 

k>0 (4 k + 3) (2k + 1) { ) 4 

2k". 
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Translation transforming the zeta 2 series, we get 

2 

(n ■+ k) 

r - 

k,n 2 

(n + 2 k) 

2 
(2 n + 5 k) (n + k) 

g - 1 + 

k,h 2 

2 (2 k - 1) (n + 2 k) 



3 3 

k (n + k) 

u - ■ 

k,n 2 2 

2 (2 k - .1) (n + 2 k) (n + 2 k + 1) 



Here we must use n « and sum from k - 1 so that r 

k,n 
is initially that of zeta 2: 

3 
21 k - 8 k 

zeta 2 - R ; 

8 (2 k - 1) . k>l 2 

8 (2 k - 1} (2 k + 1) 



3 
21 k - 8 ■ k 

R 

8 k>l 3 

8 (2 k + 1) 



S 
1 '(21 k ■•■ 13) k! 



8 j±~d 3 

k>0 (2 k ■ + 1) { 



providing B bits/term. 
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Translation transforming the 4/pi formula, 

2 . 
(2 n + 2 k - 1) 

r. ■ " ■ 

k,n 2 

(2 n + 4. k + 2) 



2 
4 n + 10 k + 7 (2 n +2 k - 1) 
s - 1 + _ 

k,n 8 (k + 1) 2 

(2 n + 4 k + 2) 



■2 2 

(2 k + 3) (2 n + 2 k - 1) 

u o 

k »n 2 2 

8 . (k + 1) (2 n + 4 k + 2) (2 n + 4 k + 4) 



Noting that s . reduces to 
k,2 





42 k + 47 




32 k + 32 


we get 






4 

— . 1 . 

pi. 



3 
1 42 k + 5 (2 k + 1) 

- + R 

4. . 2848 k>l 3 

512 (k + 1) 



2 k - 1 3 
(42 k + 5) { ) 

. - + X 

4 L~i 12 k ' 

. k>l 2 

which, like its predecessor, is also a series of exact binary 
fractions, its kth term being an integer of about Gk bits 
shifted about 12k bits to the right of the binary point, for 
a net gain of S bits/term. Jack Holloway (of the H. I. T. 
Artificial Intelligence Lab) points out that this enables us to 
compute the millionth digit of 1/p-i without computing the 
first ha 1 f mi I I ion dig i ts. 
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Translation transforming the pi formula whose ratio approached 
2/27 (where we had. set out to beat Euler's transform), 
we have 

(2 n + 3 k + 1) (2 n + 3 k + 2) 

r » • 

K.n (2 n + 5 k + 2) (2 n + 5 k + 3) 



(2 n + 3 k + 1) (2 n + 3 k + 2) (2 n + 7 k + 5) 
s - 1 _ 

k » n 2 (2 n + 5k +2) (2 n + 5 k + 3) (2 n + 5 k + 4) 
- - (k + 1) . 2k + 1 2n + 3k + 1 2n + 3k + 2 2n + 3k + 3 

U - _ ; 

k,n 2n + 5k + 2 2n + 5k + 3 2n + 5k + 4 2n + 5k + 5' 2n + 5k + 6 

which promises a term ratio approaching 54/3125. But this 
rapid convergence is marred by .the fact that s has an 

k,8 

unfactorable cubic numerator. Fortunately, s does not, 

k.l 

so we can again pull the trick of transforming the series with 
its first term borrowed. Thus 

P' , dlk + 13) (17k. + 12) . 1 2k + 1 3k + 3" 3k + 4 3k + 5 
— - 1 R 

4 ^29 k>0 5 5k + 7 5k + 8 5k + 9 5k + 11 

or 



E(Uk + 13) (17k + 12) 

5k + 6 k 
k>0 3 (2k + 1) (k + 1) (2k + 3) { ) (-2) 

2k + 3 
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Translation transforming the special cases of the hypergeometric 
function near the beginning of this paper., we get a few familiar 
identities and a few weirdos. Recalling 



n .+ a n + b 
F(a,b;c;z) - R — _ z 

n>0 n- + c n + 1 
For F(a,l;c;z) we will have 

• n + k + a 



k, n 



n + 2 k + 



Then 



k,n 



1 - 



n + k +. a z 

n + 2 k + c z-1 



(z-2)k+(c-a)z-c k + a k + c - "a z 
F(a,l;c;z) - R ^ 

c (z - 1) k>3 2k + c + 1 2k + c + 2 z - 1 

Taking the resultant with respect to k to find simplifying values of 
a and c, we find that the degree of the term ratio reduces in the cases 

+h" o a -~ 1 * 2 a * 2 3 + X * and 2 a + 2 » but oni W '" 

the case c - 2 a is the result expressible as an ordinary hypergeometric 

f unct i on: - 



z - 2 . . z . 

F(a,l;a+l/2; ). 

2 z - 2 . 4z-4 



F<a,l;2a;z) 

In the other cases, the only way to assure the n+1 in the denominator 
of the term ratio requires particular values of z: 
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2a-4 

F(a,ll 2a-2; : 

a-3 



2a-3 

F<a,l; 2a-l; ) 

a-2 



2a-l 

F(a,l; 2a+l; ) 

. a 



2a 
Fta.l; 2>va+2; — ) 
a+1 



2 
(a - 2)" 
F(a-2,2; a-1/2; ) 

(a-3) (a-1)' 

2 
(a-1) 

2 
(a - 3/2) 
F(a-1,2; a+1/2; — ) 

(a-2) (a-1) 

(a - 1) (2 a - 1) 

2 
3 (a-1/2) 
F<a,'2j a+-; ) 

2 a (a-1) 

(1 - a) (2 a - 1) 

2 

3 a 
F(a,2j a+-j '-) 

2 (a-1) (a+1) 

(1 - a) (a + 1) 
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It is a manipulative pain in the neck to confirm a fact which 
is obvious when you think about it: 

Translation transforming an s-r pair and its' 
dual produces equivalent formulas. 

This is because the orphans split off by the l's are made by 
the s's in the dual, and vice versa. 



For some reason, discoveries rarely proceed along the most direct 
line. It is only afterwards that one sees how he could have found 
his result much more easily. An example is the next three zeta 3 
formulas, which are derived in seemingly the reverse of logical order 
This presentation has two advantages over Gaussian sani t i zat ion, 
in that besides providing a datum on human problem solving, 
it shows off the power of the splitting function- 
addition formula in two unusual applications. 

3 

n 

r m . 

0,n 3 

(n + 1) 

Employing the heuristics, we try 

3 
(n + b) n ■ 
s ■ 1 + , 

8,n 3 

. a (n + 1) 

a must then be 2 to reduce the degree of u, while b = 2 zeros the 
leading coefficient in an ugly quadratic factor of u's 
numerator. The result is 



3 

(n + 1) {2 n + SJ 

r m : t ■ 

l,n 3 

(n + .3) (2 n + 3) 

Only slightly discouraged by the fact that the ratio is now 
quartic, we should obviously try 

3 
(n + b) <n + 1) 
s - 1. + ~ . 

■ l.n 3 

a (n + 3) 
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But an exhaustive application of the FACTOR and RESULTANT functions 
of MACSYflA fails to turn up any promising leads. Since 
zeta 3 formulas are so hard to find, we desperately try 



3 
(n + b) (n + 1) 
s - 1 + ,-„ . 

l,n 2 . 

a (n + 3) (2 n + 3) 

Again a must be 2 to reduce the degree of u, and amazingly enough, 
b » 4 gives 



2 
(n + 2) (n + 3) 

r - 

2,n 2 

(n + 4) (n + 5) 



a simplification to something more tike the initial ratio. This suggests 
the possibility of a sequence of ratios alternately cubic and quart ic, 
generated by an alternating sequence of splitters. If we succeed, we 
can use the addition formula to combine pairs of consecutive splitters, 
to produce the compound splitter which takes us from one cubic form 
directly to the next. Anticipating thi.s, we relabel the indices so that 
the cubic ratios have consecutive k'.s, while the quart ie ones are 
distinguished by a prime {*). Thus, this last ratio (a cubic) wi I i be 
renamed r .while the previous quart ic will be called r* 

!» n 8,n 

Now we can guess that 



2 
(n + 2 k) (n + 3 k) 

■ . r » . 

k,n 2 

(n + 3 k + 1) (n + 4 k + 1) 

and we will shortly find that 

2 
tn +.5k + 2) tn + 2 k) (n + 3 k) 
s - 1 + r 

k,n 2 

(4 k + 2) (n + 3 k + 1) (n + 4 k + 1) 

gi ves 
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2 

:n + 2 k + 1) (n + 3 k + 1) (2 n + G k + 5) 



r a 

(n + 3 k + 3) (n + 4 k '+ 3) (2 n + 6 k + 



and then 



2 

(n + 5 k + 4) (n + 2 k + 1) { n + 3 k + 1) 
s « 1 + ^_ 

k t n ~ ~~~ 

(2 k + 2) (2 n + 6 k + 3) {n + 4 k + 3) 

, indeed confirms 

2 
(n + 2 k + 2) (n + 3 k + 3) 

r - : 

k+l,n 2 

(n + 3 k + 4) (n + 4 k + 5) 

Now we combine s and s' in the -addition formula and get a 
staggering expression for the compound splitting function 
involving an unfavorable polynomial of degree six in both k 
and n. Proceeding to compute u anyway, we determine that each 
of these ungainly terms will contribute a generous twelve bits: 

2 

63 k + 78 k + 13 (k + 1) (3 k + 1) (63 k + 50) 



+ 



K.l 8 (2 k + 1) (3 k + 2) 2 

132 (3 k + 2) (4 k + 3) 

4 

(3 k + 1) (k + 1) 

U - . ; 

k,l 2 2 

.16 (3k + 4) (4 k + 3) (4 k + 5) 



4 3 2 



4 



2457Bk +641Glk +62152k +26427k+4154 (k + 1) 

zeta 3 - R 



3456 (2k + 1) (3k + 1) (3k + 2) k>8 2 

.16 (4k+5) (4k+7) 

: ■ 4 ■ 3 2 4 2 

.p™, (24570k +64161k +62152k +26427k+4154) k! (2k+l) ! 

J — — —— — — — — — — — — : _ „ 

2 



2 



*>0 . 48 (2 k + 1). (3 k + 1) (3 k + 2) (4 k + 



3) ! 
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This combination of cumbersome expression and rapid convergence 

is reminiscent of too many applications of the translation transform. 

Noting that the translation transform can be inverted, we ask 

what acceleration formula can be translation transformed to get 

thi s monster. 



I nverse 


Translation Transform 




s 

k, n 


<- 


s - 1 
k,n-k-l 




r 
k,n-k-l 


r 
k,n 


<- 


r 
k,n-k 





These relations follow from substituting n-k for n and n-k-1 for n 
in the translation transform, then solving for r and s 

k,n k, n 

In the case in question, s . ' is too large to exhibit, but the result 

k,n 

of the inverse translation transform is a barely tractable 



(n +. k) (n + 2 k) 



k, n 



(n + 2 k + 1) (n + 3 k + 1) 



2 3.22 32 

(n+k) (n+2-k) (2n +(20k+12)n +(S8k +83k+25) n+80k +150k +93k+19) 

s . ■ 1 + . 

k,n . 2 2 

4 (2 k + 1) (n + 2 k + 1) (n + 3 k + 1) (n + 3 k + 2) 



then s 



1 + 



k,l 



(k + 1) (40 k + 29) 



36 (3 k + 2) 



and 



(k + 1) (2 k + 1) 



fc.l 2 2 

' 9 (2 k + 3) (3.k + 2) (3 k + 4) 



PAGE 58 



so that 



zeta 3 



2 . ■ 4 
364 k ■ + 501 k + 173 (k + 1) 
_ R __ ___ 

k>8 2 2 

9 (3 k + 4) (3 k + 5) 



144 (2 k + 1) 



■E 



2 G 

(364 k + 501 k .+ 173) k! 



k>0 36 (2 k + 1) (3 k + 2) ! 



which is somewhat nicer and still provides nearly 3 digits/term. 
Reapplying the inverse transform, 



n (n + k) 



k,n 2 

(n + k + 1) (n + 2 k + 1) 



(k + 1) (n + 3 k + 2) n + 3 k + 1 

i » + 

k.n 2 2 (2 k + 1) 

4 (n + 2 k + 1) 



so that 



zeta 3 



30 k + 19 



(k + 1) 



R 



IB (k + 1) (2k+ 1) k>0 2 

4. (2 k + 3) 



_^ " 38 k - 11 

*-** 3 2 k 2 

. k>l 4 (2 k - 1) k ( . ) 

k 

which, at 4 bits/ term, is similar to the very first zeta 3 formula 
we got, but twice as fast. It is, however, qui te di f ferent from 
and computationally superior to grouping terms' .of the slower series 
pairwise, which produce an unfactorable. cubic numerator. 
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In general, grouping terms pairwise tends to double the 

degree of the term ratio while doubling the digits/term, 

so if you intend to evaluate an R expressed series sequentially, 

the doubled computational effort per term will counter-vale 

the greater convergence rate. Also, it greatly simplifies machine 

computation of these series when the integers comprising the 

numerators and denominators in the R expression do not exceed 

the computer's register size. 

On the other hand, Rich Schroeppel has noted that n 
iterations of pairwise term grouping will, 

n 
for large enough n, evaluate 2 terms of a sum much more rapidly 
than sequential evaluation, assuming the Fast Fourier (n log n) 
multiplication algorithm. 

Combining terms pairwise is neatly done in R notation: 

. ... (b + r (b • ■ + r ( ... 

2k ' 2k 2k+l 2k+l 

«... (b +r b + r . r (... 

2k 2k 2k+l 2k 2k+l 

or 



bRr-(b+rb )Rrr 
k k>8 k ..2k 2k 2k+l k>8 2k 2k+l 
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Pitfalls 

Inverse translation, transforming arcsin 1/2, 

2 n + 1 2 n-- 2 k + 1 
r - 

k,n 8n + 8 2n + 4k + 3 

8 n 2 n + 1 

s - - ■ ._ 

k,n 3 2 n - 2 k.- 1 2 n + 4 k + 3 

gives 

2 k +1 4 k + 4 4 k + 4 

u - . - — - — . 

k,n 2 n - 2 k - 1 2 n + 4 k + 3 2 n + 4 k + S 

2 k + 1 4 k + 4 4 k + 4 

r - ' .- 

n,k 2 k - 2 n + 3 4 k + 2 n + 5 4 k + 2 n + 7 

Noting that this simplifies for n * 1, and that arcsin 1/2 begins 

11. 
- tl + — { 1 + ... 
2 24 

so that a is 1/48, 
8,1 

Pi 1 1 4 k + 4 4 k + 4 

6 .2 30 k>0 U + 7 ■ 4 k + 9 ' 

But this is absurd -because it says that 'pi "> 3 + 1/5, since the 
R expression is > 1, What went- wrong? Recall that in the 
derivation of the transform, the kth partial sum of the orphans 
comprising the new series differs from the original sum by 
the kth partially accelerated series 



£-*i k,n 



a ■ . R r«=a uu . . . u R r 
k,l n>l k,n 8,1 8,1 1,1 k,l n>l k,n 



n>l 

which we assumed went to zero with increasing k. 
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Now as k approaches infinity, r approaches 

k,n 

2 n + 1 

16 (n + 1) 

which is the term ratio for 

-1/2 
(1 + 1/8) 

by the binomial theorem, so 



2 n + 1 1 

R - 1 - - 

n>0 16 (n + 1) 



3/2 



R 



2 n ■+ 1 2 . 

16 n>l 16 (n + 1) 3 

Thus the R expression in our error term approaches 

3/2 
2 

16 ( i - ; ) 

. 3 

while the u's which are Its coefficient form a convergent 

infinite product! Ue can, in fact, evaluate this product as 

the discrepancy between the two sides of the false equation. 
To do this, we must find out what 



1 4 n + 4 4 n + 4 

30 n>0 4 n + 7 4 n + 9 
really equals. Writing 

4 n 4 n + 4 k 



k,n 



k,n 



4 n + 4 k + 3 4 n '+ 4 k + B 
n + k 



1 + 



k + 1 



Then 



k, n 



n + k 4k + 3 4k + 5 
k + 1 4 n + 4 ■ k + 3 4 n + 4 k + 5 



(whose product over k, we note with relief, "diverges to zero" for 
n > -1) and 
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4 k + 3 4 k + 5 

1 , k 4 k + 7 4 k + 9 

But this is just the sum 

111111 

2 3 5 7 9 11 

with its terms grouped pairwise. Using 

1/2 
2 x 



arctan 



1 - x 
1/2 . 



x 



3 5 7 9 

X X XX 

— - — - — + — 
3 5 7 9 



evaluated at 1, we get 

1 4 n + 4 4 n + 4 1 pi 

38 n>3 4 n + 7 4 n + 9 2 3/2 . 

2 

Final ly, then, 



2 k + 1 4 k + 4 4 k + 4 



1/2 

3 2 pi 



n 

k>0 2 k - 1 4 k . + 5 4 k + 7 8 
This relation can be verified by multiplying by 



sin p t /4 
pi/4 



n 



4k + 3 4k + 5 



k>B 4 k + 4 4 k + 4 



to leave 



n 



2 k + 1 4 k + 3 



k>0 2 k - 1 4 k + 7 
whose kth partial product telescopes to 

2 k + 1 3 
- 1 4 k + 7 
which approaches -3/2. More generally, ue have the 
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Corrected Series Transform 



R r , " s Ru - + n u ■ R ■ r 

n > m J» n k,m k>j k,m k>j k.m n>m n 

where r is the limit of r as k grows infinite, 

n k,n 

In our particular case, II u converges for all m, giving 

k,m 

2k + l 4 k + 4 4k + 4 

n 

k>j 2k-2m + l 4 k + 2 m + 3 4 k + 2 m + 5 

as a sum of two sums divided by a third sum. I wish 1 had time 
to further investigate this connection between sums and products, 

Hoping for another product formula by dualizing the transform 
that started all this, we are very unpleasantly surprised. 



2n + l 4n + 4 .4 n + 4 

r - __ 

k t n 2 n - 2 k + 3 4 n + 2 k + 5 4 n + 2 k + 7 



3 2 n - 2 k + 1 4n + 2k + 3 

s ■ - - . 

k.n 8 k 2 k + 1 



2k + 32k-2n-l 

k.n 8k 2 k + 4 n + S- ■ ■ 

The product of the u "diverges to 8", leading to a contradiction 

k t n 

of the fact that the two series differ by a finite product. 

The problem here is very serious: each individual split is invalid! 
The quadratic growth of s is too fast for the input series, 

k,n 

which converges like a reciprocal quadratic for. all k. 

It was most unnerving and enlightening of Eugene Sal am in 
of the Draper Laboratories to. point out that any series can be 
converted to any other series in a single split, since s can 

n+1 
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always be chosen to satisfy the recurrence 

1 - s + r s = anything . 
n n . n+1 

Since each of our series transforms involves applying an infinite 
number of splits, how can we be sure tie aren't perpetrating the 
wildest absurdities? Our salvation lies in zippers. To see 
what actual ly happens to our series, imagine a zipper 
traveling along, splitting it in two. Right behind it, another 
zipper is zipping it back together, but offset so that one tooth 
(the orphan) is left hanging off the front and there i s a 
corresponding extra tooth on one side between the zippers. 
If the series were finite, this extra tooth would become the 
orphan on the right when the zipper reached the end. In symbols, 
when the zippers reach the pth term, 



e-i 



La - a s + X^ 
n . Z-r 



.1 - s ) a + s a 

n n n+1 n+1 



n>0 n=0 



+ (1 - s ) a 
P P 



+ Ev 



n>p 



If the origi na I series converges, then the last 
sum will vanish as p grows, but this i sn' t necesari I y true of 
that extra tooth, (1 - s } a . Thus, besides worrying about 

p p 

the product over k of the u , we must make sure that the 

k, n" 

limit as p grows large of (1 - s ) a approaches something 

kiP k,p 

(preferably 0} which we can some over k. At last, then, the 
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Really Truly Correct Series Transfon 



R r - s R u 

n>m j » n 



r lj +■ IT u Rr+ V^ t 

k,m k>j k,m k>j k,m n>m n / j I 



*>) 



where t - Lim (1 - s ) a as n approaches infinity, 
k n k,n k,n 

r - Lim r as k approaches infinity, and, of course, 
n k k,n 

u - 1 - s + s r. Ordinarily, s would be chosen 
k,n k,n k,n+l k,n K,n 

so that the product and the t were 0. 

k 
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Summing the Reciprocal Fibonacci Numbers 

By an involved process, I found 

k 
(-) f 
n 
r « 

k,n f 

n+2k+l 

k 
f + (-) f 
n+Ak+2 n+2k+l 

s - 

k,n (2 f + f ) f 

2k 2k+l n+2k+l 

n+k+1 
{-) f f 

2k+l 2k+2 

u - ' 

k,n (2 f + f .) f f 

2k 2k+l n+2k+l . n+2k+2 

where f , f , f , .... are the fibonacci numbers 0, 1, 1, 

8 12 



so that 



E; 



k 

f . + (-} f k 

4k+3 2k+2 (-) 

„ R ____ 

f f k>8 2 f + f 

n>l n 2k+l 2k+2 2k+2 2k+3 



k(k-l) 



k 2 

(f . + (-) f ) (-) 
4k+3 ■ 2k+2 



- Z-T-7 



9 g . . . g 
k>3 2k+l 2k-2 13 2k-l 



■ n 'In 
where g - 2 f + f - 6 + (- -) , the nth Lucas number, 
n n-1 n 6 . 

This is an extremely rapidly convergent series, yielding about 

2 
k /5 digits for k terms. 
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Similar manipulations yield 



E 

n>l 



1 

sinh nx 



_ sinh (2 k + 2) x - sinh (4 k + 3) x 

£—< sinh(2k+l}x sinh(2k+2)x (2-2cosh x) (2-2cosh 3x) . . . (2-2cosh (2k+l) x) 

k>e 

flore genera I ly , i f 

n n 

' f m x + c y 
■ n 



so that 



f -{x + y)f~xyf 
n+1 n n-1 



then we have 



k k . 

x y f 
n 

r ■ 

k,n f 

n+2k+l 

k+1 
f - (x y) f 

n+4k+2 n+2k+l 

s - . 

k,n k k+1 k+1 k 

(x - y } tx - y ) f 

n+2k+l . 

n+k 2k+l 2k+l . 2k+2 2k+2 
c (x y) (x - y ) (x - y ) 

u = -? — ■'■ " 

k,n k k+1 k+1 k 

(x - y ) (x - y ) f f 

n+2k+l n+2k+2 
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so that 



E 



n>m n 



k+l. : 
f - Ixy) f m +k 2k+l 2k+l 2k+2 2k+2 

m+4k+2 m+2k+l c (xy) {x - y ) (x - y ) 

R _ 

(x - 1) (1 - y) f f k>3 k+1 k+2 k+2 k+1 . 

ffl m+1 (x - y ) (x - y ) f f 

m+2k+2 m+2k+3 

again with the. convergence rate of a theta function. 

The author is perplexed by the fact that terms in the right side 

j j+i j+i j 

of this equation blow up uhen x = y or x = y for some 
nonnegative integer j, while the left side does not. In fact, it 
is hard to see how equality could ho id when x is fixed at 2, for 
example, and y is nearly 1. Numerical experiments, however, 
indicate that the equation indeed holds, the large term being 
rapidly eroded by later ones. This indicates that when such an 
offending j exists and is positive, the above expression, with 

R Instead of R , must be 8. 
k>j k>0 

n 
If we really do wish to sum l/(x - 1), (a so called Lambert series 
when x - 1/z), we must resort to sneaky tricks: 



1 1 

1/4 n n 1/4 n n 
n>l x. - 1 n>l (x ) -■(-!) (-x ) - (-1) 



E — " Z 



(Ixl > 1) . 

This holds because, on the right, term 2 m + 1 cancels 

term 4 m + 2, leaving only the terms where n » .4 m. 

These latter two sums avoid the division by zero when fed to the 

transform. 
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In other cases ue can avoid the problem by using the identity 



n 



c y k+1 
{ ) 



n 



n 

c y k 
( ) 

n k 

x - k ' y n 
_ c ( j u 

n n k+1 n . n 

x-cy x x-cy 



x 



repeatedly to establish the Kummer transfor 



m 



1 



n n 

n>0 x - c y 



1 k 1 y k 
k>8 (-} - - ( ) 

c x c X 



(Ixl > lyl) . 
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One Last Zeta 3 Formula 

Knopp, Chapter VIII exercise 128a, points out a relation 
between zeta 3 and the series with ratio 

3 
n 



3 

(n + 3) 



Tasting blood, we propose 

3 
n 

r - 

k,n 3 

(n + 2- k + 1) ■ 

as a candidate ratio in yet another zeta 3 transform. 
After guessing the form 

4 3 2 
an +bn + c n + d n + e 
s - 1 + 

k,n 3 

(n + 2 k + 1) 

we find (with just a little help from MACSYriA) 

1 2 k ' + 1 
a ■ - , b « 

6 k + 2 3k + 1 



3 
(2 k + 1). 

C ■ : ._ 

2 (3 k + 1). (3k + 2) 

3 
<k + 1) (2 k + 1). 



2. (3 k + 1) (3 k + 2) 



2 3 

(k + 1] (2k + 1). 

e - -____ — : 

3 (3 k + 1} (3 k + 2) 
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Despite these composite coefficients, the author has so .far 
failed to find any nicely factorable expression or sum of 
two expressions for. this splitting function. Putt i n.g. thi s 
function in some intelligible form could shed much light on how and 
when to look for such things. As it stands, the gross disparity 
in elegance between the ratio and its spli tter seems to suggest. 
that we should seek a new formulation of the transform process. 
In any case., taking partial fractions with respect to n. 



3 2 2 

n + 4k + 1 (2k + 1) (3n + (15k + 9) n + 28k + 25k + 8) 

s . + 

k,n 2 (3k + 1) 3 

6 (n + 2 k + 1) (3 k + 1) (3 k + 2) 

Then 

2 . 2 
56 k + 80 k + 29 (k + 1) . 
zeta 3 - R 

3 k>8 3 (3 k + 4) (3 k + 5)' 

24 (2 k + 1) (k + 1) 

In the previous exercise, Knopp indicates that 

1 2 . 

E— — 18 - pi 
3 3 
n>l n (n + 1) 

This is 1/8 (- the first term) of what we would get. by substituting 
k + 1/2 for- k in the last R expression for zeta 3, since 

3 
n 
r * . 

1/2, n 3 

(n + 2) 



Euler's Constant 

In the American Mathematical Monthly, vol. 7B #3, 
HarB3 p273, 1. Gerst and H. F. Sandham give 



112 2 2 2 3 3 4 
gamma ■ --- + --- + _-_ + __ t#t _ + 

2 3 4. 5 6 7.8'"' 15 16 



the numerator increasing by one each time the denominator reaches 
a power of 2. Ue can write this as the double sum 
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1111111 

- . + + --- + -... 

2 .3 4 5 6 7 8 

11111 
+ + + -... 

4 5 G 7 8 



1 1 1 
+ + „ . ... 

8 9 IB 



n 
which is the series for - log 2 started at the 2 th term, 
summed for n.> 8. Recall the application of the translation 
transform to 



, 2 n - 1 

r - 

k,n 2 n + 4 k + 1. 

(the "embarrassing" formula) to get pi/4, since for k » 8 this is 
the term ratio for arctan 1. Replacing n by n + 1/2, and redefining 

n 

K,n . n + 2 k + 1 

we have, for k - 8, the ratio we want. Replacing n by n + 1/2 
in the corresponding splitting function and u, 

2 n + B k + 3 

s - 

k,n 4 (n + 2 k + 1) 



(2 k + 1) (2 k + 2) 

U m . 

k,n 4 (n + 2 k + 1) (n + 2 k + 2) 



which says 
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n>m 



n 

<-> ? m + 6 k + 3 (2 k + 1) (2 k + 2) 

R _ _ ^ 

A m (m +1) k>8 4 (2 k + m + 2) (2 k + m + 3) 



E(2 m + 6 k + 3) <2 k) ! " 

k+1 
k>0 4 ■ (a + 2 k + 1) ! 

(for m even). Thus, to get Eulers constant, we must 

sum the R expression with m = 2 , for j > 0. 

j+1 

E2 + 6 k + 3 
\ A 

j>B k>3 2 ■ + 2 k + 1 j+2k+2 j 

C ) 2 (2-1) 

2 k 
This prospect grows 

cheerful when we notice that after a few j, the R expressions wi I I 
converge with great rapi di ty, returning 2 j + 2 bits/term.- This 
is because in the term ratio, there are two m's-in the denominator 
and none in the numerator. Thus b bits of gamma will require about 
(fo-j)/(2j+2) terms of the jth series. This. says that the total 
number of terms of all series will be about 

b b-1 b-2 B b 

■- + + -—...— ^ - (In b - 1) 

2 4 6 2b 2 
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Zeta 2 Revi si ted 



n (n + k) 

r » 

k,n 2 

(n + k + 1) 



n + 2 k + 1 . 

s • ■ 

k t n k + 1 

gives 

(k + 1) (n + k) 

u - 

k,n 2 

(n + k + 1) 

Now s « 2 and 

2 
fk + 1) 

u = 

k,l . 2 

(k + 2) 

which is a strange way to prove zeta 2 = 2 eta 2 . 

Translation transforming, 

18 k - 3 k 
zeta 2 » R 

. 4 k (2 k - 1) k>l 2 (2k + 1) ". 

This series has similar convergence to our first (Markov) on< 

3 . k 

zeta 2 = — R. -: 

2 k k>l .2 (2 k + 1) 

except that the former one alternates. 
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The zeta 2 term ratio 

-.2 
n 

2 
(n + 1) 

i s a speci a I case of 

n ■+ p n + q ■ 

n + r n + s 
Choosing the splitting function a n + b will produce a new ratio of 
n 4- p n + q . Q(a,b;n+1) 

n + r + 1 n + s + 1 Q(a,b;n) 

Where Q is a polynomial quadratic in n and linear in a and b. 
Now a and b can be chosen to make Q divisible by any two 
of the quant i ties n + p,n + q,n + r,or n '+ s , or 
they can be chosen to zero one or both leading coefficients 
of Q. Thus for our new ratio we may have any one of the 
quant i ti tes 

n + p n + q n + p . n + q + 1. 

■ — — _________ ^ ■ ___ . _ __ __ ^ 

n + r + 1 n + s + 1 n + r + 1 n + s + 1 

n + p n + q n + p n + q 

—————---- — -- , _____---_ __________ t . 

n + r + 1 n + s n + r + 1 n + s-1- 

n + pn + q n + p + ln + q + 1 

n + r n+s n + r + 1- n+.s + l 

or any of the above with p and q swapped or r and s swapped. 
Also, from the addition formula, we learn that p - 1 and 

n 

q - - 1/r combine to produce 8, the identity split, 
n n 

and thus the sp I \ tting function 

n + r - 1 : n + s- - 1 

n + p-ln + q-1 

will perform a reverse translation (the opposite of s =1) 

n 
producing 
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n + p-1 n + q-1 

n + r-1 n + s-1 

By combining the above operations, we see that it is possible 
to independently change any of the four terms of our original 
term ratio by any integer. ■ Thus, any series whose ratio is of 
this form, with p, q, r , and s all integers, must evaluate to 



c pi + d 

with c and d rational, since the splitting process only adds and 
multiplies by rational quantities.. Among these rational quant it 
is 0/0, which figures into those splits which would result 
in cancellation to a linear term ratio.- For example, to get 

n + p n + q 

n + r + 1 n + s-1 

we need the splitting function 

s-l-r n + s-1 



i es 



s-l-p s - 1 - q 



Thus, i f both denominator terms exceed both numerator terms, 
there is no way to pass through any ratio where a 
denominator term is as large as a numerator term. In this 



way we are prevented from proving that pi is rational, 
since if is easy to show that any series with ratio 

n + p 



n + r 



telescopes to a rational sum when with r is rational and 
r - p is an integer. As a corollary, we observe that we- can 
achieve such a cancellation if a denominator term initially exceeds 
one »n the numerator, thus any series with ratio 

n + p n .+ q 



n + r n + s 

sums to a rational if p is rational, r 
q - s is a nonnegative integer. 



- p i s an integer, and 



A- particularly neat special case of the above 
i s the transform; ' 
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k,n (n + c) (n + k) 



n +. c - 1 



k,n k + c - 1 



k,n (k + c - 1) ( n + kl 
which says 

2 



R 



m + c - 1 



n>m (n + c) (n + j) j + c - 1 k>j (k + c) (k + m) 



i . e. 



(j + C - 1) R 

n>m (n + c) (n + j) 

is a symmetric function in c, j, and m. 

With ail of these relations, there must.be a neat irrationality 

2 
of pi proof in here somewhere. 
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For those readers who have found the treatment so far a bit on the 
Heaviside, warning: we are approaching a singular i ty~we will now 
attempt to accelerate the "series" 



e E (1) 
1 



B!-l!+2!'-3!+... 



1-t 
gotten from iteratively integrating by parts e dt/t from 1 

to inf ini t.y. 

Ue have 



l.n 



- n . 



This is a special case of 



k,n 



n (n + c) (n + d + 1) 
k 

(n + a + k) (n + d ) 

k 



with d » c ■» a and z ■ -1. 

1 

- 1 



k,n z (n + d ) 

k 



spi i ts thi s to yield 



(1 + z td + k"+ a - c)) n ' + (k + a) (1 .+ z d ) 



k, n 



(n + d ) (n + a + k) z 
k 



which formally satisfies the recurrences if 

(k + a) (1 + z d ) 
k 



k+1 1 + z (d + k + a . - c) 
k 



Thi s says that 
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where 



- 1 
R r -■-- — „ R v 

n > m J» n z (m + d ) k>j k,m 

J 



1 + z (d + k + a - c) 

k 
v ■ «• 

^.n (n + a + k) z 



1 + zJ2 k+2a-c-l) (k +■ a - 1) (k + a - c - 1) 

tn + a + k) z "(n ■+ a + k - 1) (n + a + k) v 



k-1 



(using the recurrence relation for d to eliminate it). 

k 

Expressed in standard notation, 



V^ a •'-" a V^b 
Z-*i j»n j,m / j k,m 

n>m k> j 



a b 
k '"+l k+l,n 
Wlt h r « and v . — 

k » n a k,n b 

k »" k,n 



so that b 

k+l,m 



I '+ z (2 k + 2 a - c - . 1). , . k + a - 1- k + a - c - 1 
b _ ^ 

(m + a + k) z k,m ■ m + a'+ k m + a + k - 1 k-l,m 



with b - s and b - s • v . 
j» m i*" 1 j+l.m j,m j,m 

In the special case of 0! - 1! + 2! - 

we have j-m-1, c-a-d - arbitrary! and z - -1. 

■ 1 
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b 
k 



Choosing c - 1, we have 



n>0 k>l 

with b" - 1/2, b - 1/6, and 
1 2 

2 k - 1 . k - 1- k 

■ k+1 k + 2 k k + 1 k + 2 k-l' 



For reasons unknown to the author, this last series seems to 
oscillate with ever decreasing frequency and amplitude about 
the "correct" answer .596347..., gotten from the integral. 
The convergence is slow: summing through b gives 

.76 

.593299... and through b gives .5935.... . 

IBS 

Larger values of c tend "to reduce the frequency of early 
oscillations and the amplitude of later ones, slightly improving 
the long term convergence. 
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Numerical Experiments 

The reason that we choose s approximately - 1/(1 - r ) 

k,n k.n 

and not. exact I y is that r usual ly becomes an explosively 

k, n 

complicated expression after a few k. But there is nothing 
preventing us from numerically grinding out a few orphans 
to see how fast the series would go, could we find its 
formula. Ue try good old 

pi - 4 (1 - 1/3 +1/5 - 1/7 + ... 

and get the following data. (The orphans are consecutively 
numbered, with the cumulative sum on the next line. The correct 
portion of each sum is underlined. Sums and orphans have all been 

69 

scaled up by a factor of 10' and truncated to integers.) 

0. 30008800000008.0000030800003000000000803000000000030300008000.0000000300 

3000030000300000000080000888000000000300000000000833080030200000000030 

1 13888S8S8S88888S388888838S88SS8S8838388888888888838S83888838888888S83 
313888888888888888888388SS8888888S888888S8888838S88SS88S88S8888888SS88 

2 .26721,G184977729S94S43822364477484376233S84993217062G8391949G714G490 
3141S610S07386G6184834327112S336637326Si884738S210S9S14838B838S683S378 

3 313S2227689441396G37992SSSS183206S0334026S7822222369422124478S286 
31415924029663SS62623096S13S092182BS330221876S399281743B9238S100320G64 

4 24939GG112473911216843G792G1111891GGG3423G2G22170588G5571729G97 
3141S926523G296687362208678946010716449138S39963SS433636S1170672SS0361 

5 ■ 123173417G3487814445G59S377171482811093G41308G52G1081 11308129 
3141532G535947810499708G823402GB3893G28G213B1857195784717727878385S49B 

G - 7888184G397176797723BB514G93761B884288812B244B73BG3B9B4831 

3141B92GB3B869088GB3311BB5B42B3740421B124B2422143075844731S43147903G59 

7 . 2891758706224834123GB323381G34488032393434BG48B4494309G724 
3141B926B3B89792G2403737B388377343S2B32879G4224G70G81SB3800830970003S3 

8 32373B4G77832473G0G9B814BGB1G41323G374B064B733GG47734S 
3141S926S3S.89793S44822843I716247013211433S6S8GS99423925444G674G3477421 

3 ■■ - 307978293062G8777178826B18G618483534771G4222B7B04GS79B 

3141B92GB3B837932368445B010893G929B328782B00G811B88915380244171300SG3G 

10 17G74088051379GG711BS970B978848491072884G1712754450B 
3141B32G53B8979323SG119BS9140748962444379BG04G964380226G870B8840553141 

11 - 24127G143030B42S8439G7B91G302B9240G113117B0582349G4 
3141B92GB3B89733238370G827710442B33B944119S88333845G1GB43BS83782318177 

12 " 3170B7998G4G3348204737373BG239143G979GG10114571337 
31415326B3B837332384G2388S7890893684148SS78619SS237B22B23249389G889434 

13 218172S7880486S1BBG39121231284G130S427248B8G867 
3141S926S3S897932384G241838S29S8S732S304139S3168S433696G3832114S4763S1 

14 _ - 25S4314783497B43GB082730137323238114G4818370 

3141B92GB3B897932384G24103G24B3720233328783B234124849S4S77189G88G57981 

15 23298B22488GS9.92777812938778G4B3S8G32774B9B21521 

3141B92GB3B897932384G2G43342G78G07132GBG484314290249B314453871396S9B02 
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1B ^ 3011768197SG15258825723637134583G39G455S1667 

314159265358979323846264227279620311122174314031395663632^826785241169 
17 ' A 842452315262727517163061437334656142244305 

31415926S35897932384626433736336S192G48447B6571820181339636S2927485174 
18 01 , , ___ ~ 47883177665521S71424S5649B86672905B605936 

3141532653S837S32384626433735S1778748813S43B357S7461G370368S387GS732S8 
19 , , rM ^^ ftM „^ 968S65976B934822448413828214653S38G58783G02 

3141592653589793238462643383288438509753771534178439851755792535662878 
2 ® ,™ „_„_ " 2212139288G5692447844914358845303680 

31415926S3S8979323S4626433S32S843828S539S4266847SB5230584B933G93262190 
2l 01/1roonn , roo _ _ " 1808933343637333331178873274388527 

. 3141592G535837S32384G2G433832S843S28G730S4331SS48671G756G20G0415373GG3 

22 ^ ___ „ ' 337363528435513888425916164736 
3141592G535897932384G2G433832884382S673118GG823G91671894784SG332138459 

23 - 935556387196828079374143488869782438074 
31415926535897932384626433832795827384229898542891 

24 „ 26873347497834B433S915317272G7086 

31415926535897932384626433832795027386847703232595 

25 „ w « ,-«„„,.„_ 776379841139153325G33G2897853SG 

31415926535897932384626433832795027306925341276789 

2S _ ,,_ - 2523054G485251'21915223G58804 

31415926535897932384626433832795827336925316376162 
27 ' , 2261723SG8383G85340815055G3540S33315 

31415926535897932384626433832795829568648377859848 ■ 

2 * , _ -1491358613247912513790835335615 

3141592653S897932334626433832735B29S6S634'0574736'S6 

29 ,.„ 38633836935074657867266587 

3141592G535897932384G2G433832795B295G8G34B57779987 

3B „/,r m rr^ n , " 72G73 ^733818B031 715841434344198353 

314159265358979323846264338327950288419 

31 ow,m W nr B ^ " ^386665699367779252120813.98108 

31415926535897932384626433832795328841 

32 ,,_„™ • ■ 9816226177585303871383133436275 

314153265358373323846264338327350288413 

33 314159265358973323846264338327353288419 2 ^389186933100428338 
34 31415S2G5358979323846264338327958288419 222 SS32222G73517598341 

7O0CO1 COA9C i 7 

. 314159255358979323846264338327950288419 

36 01 , , ^^_„ - 1246586845573268956321132851058 

31415926535897332384626423832795028841971633 

37 - 24436702417395191 

■31415S2G535897932384G2G433832735028841371693- <^'^i AttbiiW 
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Ue note several things: 

1) The first. six terms are positive, after which the signs 
are apparently random, 

2) The series converges faster and faster until the 
first negative term, which provides only one 
additional correct digit. 

3) After this, the terms grow erratically smaller until 
term 15, which has nearly 1033 times the magnitude of term 
14 and yet provides four new correct digits.. 

4) Terms 18. 22. 26, 29, 31, 33, 34, 35, and 37 have 
the wrong sign, i.e. they actually make the sum worse. 
Term 31 actually undoes a good digit. 

5) . Term 3S is seventeen orders of magnitude larger than 

term 35 and provides five new digits, while the sizable term 27 
overcompensates and contributes less than one. 

6) The general trend after term 6 seems to be ever slower 
convergence. 

Little wonder that we couldn't find a formula — I dare 
you to prove i t converges at all! 

I am indebted to Richard P. Howell of the HIT RLE PDP-1 
for. confirming that these results are not merely the PDP-18 
bignum routines gone haywire. 

Now in the early 1958's, Oaniel Shanks did a- Ph. D. thesis 
(Journal of Mathematics and Physics,. 1955, first article) 
involving primarily numeric methods. of sequence extrapolation. 

Among them was a transform which he denoted e , 

1 

which, in effect, replaces each but the first term with 

the limit of the exponential sequence defined by that term's 

predecessor, itself, and its successor— the first member of the 

kth sequence produced by this process becoming the kth member 

of the transformed sequence. On page 5 Shanks gives a 

numerical example for which he too chose the 4 arctan 1 

series (actually the sequence of partial 

sums). He computed only six transform terms to eight decimal 

places, since he was technologically constrained to electro.- 

mechanica! computation. Rich Schroeppel made the remarkable 

observation that if Shanks had started his sequence with 

a zeroth partial sum of 8, his transformed sequence would have. 

been identical to the screwy sequence of partial sums 

that we got from splitting by l/(l-r), and indeed, extending 

the computation of his example (without the initial 3) 

exposes the. same erratic convergence. Although "it is doubtful 

that he anticipated this, Shanks provides some valuable insight 

into this phenomenon with his section on spectra of sequences. 
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Miscellaneous Accelerations 

In. exercise 129c of chapter VIM, Knopp gives a relation 
completely equivalent to one step of a transform for eta 2 and 
Catalan's constant, but he seems to miss the significance of it, 
presumably because he was more interested in identifying series 
than in evaluating them. Rephrased in terms of splitting functions, 

2 
n 

r « . 

k,n 2 

' (n + 2 k) 

i k 2n + Gk + l 

k,n 2 2 .2 

(n + 2 k) 

3 
k (2 k + 1) 

u - . _ 

k f n 2 2 : 

2 (n + 2 k) (n + 2 k + 1) 

Inverse translation transforming the "embarrassing" pi/4 
formula: 

n - k 

r - 

. k,n n+3k+l 



1 (4 k>l) (4 k + 5) (4 k + 1) (2 k + 1) 4 k + 1 

S m - + -i ;_ + 

k » n 2 16 (n + 3 k + 1) 8 (n + 3 k + 2) 16 (n - k - 1) 



u 



1 k +- 1 2 k + 1 4 k + 1 4k + 3 

k » n 2 n-k-1 n + 3 k + 1 n + 3 k + 2 n + 3 k + 3 

so if n starts at 1/2 and k starts at 0, . 

P' llk+8 8k +1 4k+14k+3 
— , R _ _ : _ ___ 

4 15 k>B 3 2 k + 1 6. k + 7 6 k + 11 



£ 



k 
{11 k + 8) 16 k! (3 k) ! (4 k) ! 



3 (6 k + 5) (2 k) ! (6 k + 1) ! 

k>B 

which, with a term ratio approaching 16/27, - has the least - 
finite convergence rate of any (rational) pi series I know. 
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Even worse is the bits/term one gets by purposely choosing 
sp I i tting functions which don't make u tend to zero with n: ' 

"n + k n-k + 2 
r « ■ '. 

k,n n + k + 2n~k + l 

2 n + k 
s .1 

k,n 4 k + 3 n-.k + l '. 

4 k - 1 n + k n - k 

u - ■ 

k.n 4k + 3n + k + 2 n-k + 1 

Choosing n - -1/2, k - (i.e. starting the arctan 1 series one 
term too ear I y) , 

pi 4 k + 5 ■ 2k + l U-'l 

4 4 k + 3 k>8 2k + 3 4k + 3 

k 
(4 k + 5) (-1) 

Zj (2 k +1) (4 k - 1) (4 k + 3) 

k>0 

which is no better than i f we had taken terms of the original 
series pairwise. Grouping pa i rwi se, however, would have 
destroyed the a I ternati'ng. sign. Instead, this series 
arises directly from arctan 1, 

2 k' + 1 , 
r - - , 

. k 2 k + 3 

with the single application of the splitting function 

2 k ' + 1 

s - .- ... 

k . 4 k - 1 ■ 



Retransforming Euler's transform of arctan 1, 



0,n 2 n + 1 



I found 
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n + k 



2 n + 4 k + G k + 3 



k t n 



k,n 



2 n + 4k+'3 2 

2 n + 4 k + 6 k + 1 



4 n + G k + 3 
2 n + 4 k + 3 



2 k - 1 



n + k 



2 n + 4 k + 14 k +'11 



k, n 



2n + 4k + 3 2n + 4k + 5 2 

2 n + 4 k + 6 k +■'! 



which is very unlike the other accelerations in this paper 
due to the quadratic occurrence of k in r and u. 

Translation transforming the foregoing, then initial i z i ncj 

n to 1 and k to 0, 



pi 8 (14 k + 33 k + 20) 2 k + 2 2 k - 1 

2 1B5 ■ k>8 . 3 B k + 11 ' S k + 13 

Noting" the similarity between this formula and our earlier 

Pi 2 (7 k + 8) . k + 1/2 k + 1 1 

4 15 k>0 k + .7/6 k + 11/B 27 



ue suspect a single application of some splitting function relates 
them. To find it we divide, the first formula by 2 and rewrite the 
second formula. so that they both have the same expression to the 
right of the R. 



2 (7 k + G) 

15 ■ 



R 

k>0 



2 k + 1 2 k + 1 k + 2Sk + 13 2 k. - 1 

3 G k + 7 .6 k + 11 k-+ 2 G k + 13 "2 k - 1 



2(7k + G) 2 k - 1 Gk + 7 i 



15 



- 1 



k+1 k>0 



2 k + 2 2 k- 1 

3 G k + 11. G k + 13 
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b b 

• k+1 k 

(Thrice using the rule -R r - — R r 

k>0 b ■■ k b k>0 k 
k 



Now we have two series of the for 



m 



b R r and c R r . 
k k>8 k k k>0 k 

Performing a modified split on the first of these, 

(b - t + t J R r - t ■ R r + (b - t ) R r . 

k k k k>0 k k k>0 k k k k>0 k 

- t + "p t R p + (b - t ) , R p 
8 k+1 k>0 k+1 k k k>0 k 



t + p f ftp + (b - t ) R r 

8 k k+1 k>0 k k k- k>8 k 



-■ t + (b - t + p t ) R r 
8 k k ■ k k+1 k>8 k' 

which we must' make equal to c R ' r with the right choice of t . 

k k>0 k k 

This means t '- 8 and 

t - b- + c ■ . 
k k k 
t « _ , 

k+1 r 

k 

In thi s case, 

2 

4 (14 k + 33 k + 28) 
b » . 

k 105 

2 (7 k + G) (2 k - 1) (6 k + 7) 

■ c - 

■ k 105 (k + 1) . 



2 k + 2 2 k - 1 
k 3 6 k + 11 S k + 13 
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Factoring the first feu numerical values of t , we can easily guess. 

k 

2 k 6k-+5Gk + 7 



k k + 1 S 



7 



which is confimed by direct substitution into the recurrence. 
Noting 

(b - t ■ + r. t ) R r 

k k k k+1 k>8 k 



t t r b . b 

k k+1 k k+1 k+i 

b (1 - — + ) R r 

8 b b b k>0 b k 

k k+1- k k 

we see that t /b is the ordinary splitting function 

k k 

for the series with ratio b r / b . 

k+1 k k 



Thus 



3 k (G k + 5) (G k + 7) 



2 <k + 1) (14 k + 33 k + 28) 



spl i ts 



pi 4 (14 k + 33 k + 2B) 



105' 



2 k + 2 2 k - 1 

k>a 3 G k.+ 11 G k + 13 



to get 



Pi.. 



2 (7 k + S) 



15 



. 2 k + 1 2 k + 1 
k>0 3 G k + 7 6 k + 11 



eliminating that ugly quadratic. Perhaps these bothersome ..■' 
polynomials which" arise so often in our faster accelerations are 
si mi I'ar I y dispensable. Perhaps there is even a way to do this 
to series whose convergence rate has been doubled several times 
by pairwise term grouping. 
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Negating both sides of the recurrence for t 

k 

. (- t ) ■- c + b 
k ■ k k 
{- t ) » 

k+1 r 

k 

.which says that if t takes b to c , then -t takes c back to b . 

k k ■ k . k k k 

Thus 

3 k S k + 5 
k 7 k + G 2 k - 1 
carries the second series back to the first. 

Another example. of alternating splitting functions: 
First a pre I i mi nary spl i t of 

n 

r m l_ 

n 2 n + 1 

by 

4 n + 4 

S m 

n 2 n + 3 

to get 

n 2 n - 3 

0,n 2n + 7 2n-5* 

Then 

n + k 2 n -. 4 k - 3 

r - • . 

k t n 2n-4k-5 2-ri + 8k + 7 

i s spl i t by 

4n+'4k+4 
p - 

k,n 2 n + 8 k + 7 . 

to get 
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n + k 2n-4k-l 
r' - - 

k,n 2n~4k-3 2n+8k+ll 

which is in. turn split by 

4 n + 10 k + 11 : 
q _ 

k,n 2 n + 8 k + 11 

to get r p and q combine in the addition formula to 
k+1 

2 n - 4 k - 3 Bk + 5 Gk + 9 
s - p + q _ 

k.n k t n k,n 2 n - 4 k ~ 5 2 n + 8 k + 7 2 n + 8 k + 9 

which gives a u containing five linear terms over five linear terms 

k,n 

which approaches -27/512. s is quartic, and for at least the small 

k, n 

positive integers, won't factor. The resulting pi formulas are dull. 
Mildly interesting is the fact that translation transforming 
yields a formula both faster and slightly simpler, but not enough 
to make i t a good pi formula. (Unfactorable cubic on the left of 
R, four linear factors, over four linear factors on the right, 
ratio approaches -27/3125.) 

Somewhat different pi formulas result if the preliminary splitting 
function s is merely. 2, which gives 
n 

n 

r « » r' 

. 0»n 2 n + 3 -l,n+l 

so that q and p ' can be combined with the addition 
k-l,n+l k,n+l 

formula to get a different splitting function. 

In general, the arguments to the addition formula my be 
interchanged by making the first argument a preliminary split 
on the first iteration. Unless one of them is identically 1 
(the translation transform), the resulting formulas will 
differ significantly, but .converge at the same rate. ' 



